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ABSTRACT 

The  extent  of  molecular  mixing  in  several  two-dimensional 
free  turbulent  shear  flows  was  measured  using  a concentration 
probe  with  a frequency  response  of  100  kHz  and  a spatial 
resolution  of  0.  1 mm.  The  flows  investigated  were  (i)  a shear 
layer  in  which  the  gases  on  either  side  of  the  layer  are  of 
unequal  density,  (ii)  a shear  layer  in  which  the  gases  on  either 
side  of  the  layer  are  of  equal  density,  and  (iii)  a wake  in  which 
the  gases  on  either  side  of  the  wake  are  of  unequal  densities. 

The  extent  of  mixing  was  measured  as  a function  of  Reynolds 
number  for  the  first  case. 

It  was  found  that  at  a critical  Reynolds  number  the  extent 
of  molecular  mixing  sharply  increased  (25%).  Power  spectral  density 
curves  of  the  concentration  time  histories  also  indicated  a marked 
increase  in  the  high  frequency  fluctuations  above  this  Reynolds 
number.  A shadowgraph  investigation  of  this  phenomenon  revealed 
that  three-dimensional  Taylor -type  vortices  whose  axes  of  rotation 
are  basically  in  the  flow  direction  exist  in  the  flow  in  addition  to 
the  two-dimensional  large  structures  previously  observed. 

These  Taylor  vortices  were  found  to  be  unstable  above  the  critical 
Reynolds  number  and  were  producing  the  increase  in  molecular 
mixing.  The  growth  and  development  of  the  two-dimensional 
large  structures  were  found  to  be  basically  unaffected  by  this 
instability.  It  is  proposed  that  the  fully  developed  turbulence 
of  shear  flows  is  maintained  by  a combination  of  the  development 
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of  the  large  structures  and  of  the  coupling  between  the  large 
structures  and  these  unstable  Taylor  vortices. 

These  data  were  also  used  to  predict  results  for  shear 
flows  in  which  diffusion -limited  chemical  reactions  have  been 
incorporated. 
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I.  INTRODUCTION 


This  experimental  investigation  had  its  prime  motivation 
from  a continuing  effort  at  the  California  Institute  of  Technology 
directed  toward  the  understanding  of  turbulent  shear  flows  in 
general  and  turbulent  and  molecular  mixing  in  such  flows  in  parti- 
cular. 

Within  this  general  framework,  the  mixing  process  in  a 
two-dimensional  turbulent  shear  layer  produced  by  two  fluids  (both 
of  equal  and  different  densities)  flowing  parallel  to  each  other  at 
different  velocities  was  studied  (see  Fig.  1).  The  aim  was  to 
determine  the  extent  of  molecular  mixing  within  the  layer  and, 
further,  to  determine  what,  if  any,  effect  Reynolds  number  has  on 
the  extent  of  the  mixing.  (The  possibility,  in  fact  probability,  of 
such  an  effect  is  evident  from  Plate  1 which  shows  shadowgraphs 
of  this  flow  at  various  Reynolds  numbers.  In  all  photographic 
plates,  the  flow  direction  is  from  left  to  right  with  the  high  speed 
side  on  top.  The  conditions  for  all  photographs  are  listed  with 
each  plate.  ) 

This  flow  geometry  was  chosen  for  several  reasons.  First 

of  all,  the  flow  is  well  documented  both  for  the  case  when  the 

1 2 

two  fluids  are  of  equal  density  ' and  for  the  case  where  the  two 
fluids  are  of  unequal  densities^’  Secondly,  this  mixing  process 
is  similar  to  the  process  by  which  two  reactants  are  mixed  in 
many  continuous -flow,  nonpremixed  chemical  reactors  or  combustors. 
However,  at  present  little  analytical  or  experimental  data  are 
available  which  can  be  used  to  predict  the  extent  of  mixing,  and 
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thus,  reaction  rates,  even  for  a shear  layer,  much  less  for  the 
often  complex  flow  dynamics  in  practical  reactors.  For  compari- 
son, the  mixing  in  a two-dimensional  wake  was  also  studied. 

A further  and  possibly  more  important  aim  of  this  investi- 
gation was  to  shed  some  light  on  the  continuing  discussion  among 
fluid  mechanists  concerning  the  fundamental  phenomenological 
laws  which  govern  free  turbulent  shear  flows.  Until  recently  it 
was  thought  that  such  flows  consisted  of  two  regions,  one  turbulent  and 
the  other  nonturbulent,  with  the  turbulent  region  being  character- 
ized by  random  three-dimensional  motions  and  the  presence  of 
vorticity  fluctuations.  Viscous  forces  were  continuously  propa- 
gating these  vorticity  fluctuations  into  the  nonturbulent  region 
along  the  interface  between  the  two  regions  (entrainment  by 
nibbling).  By  this  process  the  turbulent  region  was  continuously 
being  enlarged  at  the  expense  of  the  nonturbulent  region.  Due  to 
the  recent  discovery  of  large  coherent  structure  in  many  turbulent 
shear  flows,  a new  view  of  turbulence  is  now  emerging5'10. 

This  view  suggests  "that  with  every  shear  flow  is  associated  an 
identifiable  characteristic  structure  and  that  the  development  of  the 
flow  is  controlled  by  the  interactions  of  these  structures  with  each 
other10".  Thus,  an  understanding  of  the  development  of  these 
structures  and  their  interactions  with  each  other  should  give  some 
insight  into  the  actual  physical  processes  in  turbulent  flows,  such 
as  entrainmenl,  transport,  mixing,  noise  production,  gustiness, 
etc.  , and  should  lead  to  improved  methods  for  analyzing  and 
computing  them  . However,  ideas  concerning  the  forms  of  the 
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large  structures  in  each  flow  differ  greatly.  In  fact,  this  view  of 
turbulence  itself  is  far  from  being  generally  accepted.  For 
example,  it  has  been  stated  that  the  basically  two-dimensional 
large  structure  of  a shear  layer  as  reported  by  Brown  and  Roshko^ 
and  others  cannot  undergo  vortex  stretching  and  therefore  cannot 
form  part  of  the  fully-developed  energy  cascade  from  mean  flow 
to  smallest  eddies  which  is  the  essential  characteristic  of  turbu- 
lence. Further,  these  isolated  large  structures  when  subjected 
to  small  perturbations  would  deform  indefinitely  by  self-induction, 
thus,  leading  to  a rapid  breakdown  to  three-dimensionality.  These 
are  certainly  valid  arguments.  Therefore,  a further  aim  of  this 
research  was  to  reconcile  these  views  of  the  large  structure  in 
shear  flows  and  thus,  -to  add  to  the  basic  understanding  of  the 
nature  of  turbulence. 

With  these  goals  in  mind,  the  investigation  proceded  as 
follows.  Using  an  improved  version  of  the  Brown-Rebollo ^ concen- 
tration probe,  the  extent  of  mixing  in  the  turbulent  shear  layer 
produced  by  parallel  streams  of  high  speed  helium  and  low  speed 
nitrogen  was  determined  as  a function  of  lateral  position  (y) 
and  Reynolds  number  (Reynolds  number  based  on  AU,  vorticity 

thickness,  and  the  kinematic  viscosity  of  nitrogen  varied  from 
4 4 

0.  3 x 10  to  7 x 10  ).  The  experimental  facility  and  instrumen- 
tation used  are  discussed  in  Section  2.  The  actual  procedures 
for  data  acquisition  and  reduction  are  discussed  in  Sections  3 and 
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For  comparison,  two  other  shear  flows  were  also  studied. 


In  one  flow,  the  extent  of  mixing  in  a shear  layer  between  streams 
of  equal  density  fluids  was  measured.  For  this  flow,  nitrogen 
was  used  as  the  high  speed  gas  and  a mixture  of  argon  and 
helium  arranged  to  have  the  same  density  as  nitrogen  was  used  as 
the  low  speed  gas.  In  the  other  flow,  the  velocity  of  the  nitrogen 
was  adjusted  to  be  equal  to  the  velocity  of  the  helium,  thus 
reducing  the  shear  layer  to  a wake,  the  wake  of  the  splitter  plate 
which  initially  separated  the  two  fluids.  Thus,  the  extent  of 
mixing  in  an  inhomogeneous,  two-dimensional  wake  was  investigated. 
The  results  and  comparison  of  these  three  cases  are  in  Section  5. 

These  measurements  of  the  extent  of  mixing  in  turbulent 
shear  flows  revealed  a strong  Reynolds  number  effect  on  the 
mixing.  It  was  felt  that  this  was  due  to  an  instability  within  the 
large  structures.  In  order  to  confirm  this  idea  and  also  to  study 
the  three-dimensionality  of  these  flows,  a shadowgraph  investi- 
gation was  made.  These  results  are  presented  in  Section  6. 

Although  the  flows  studied  are  nonreacting,  the  results 
were  used  to  predict  the  reaction  rates  for  an  interesting  md 
useful  case  in  which  flows  containing  chemically  reactive  species 
can  be  successfully  studied  analytically.  The  results  for  the 
three  cases  are  in  Section  7. 
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II.  EXPERIMENTAL  FACILITY  AND  INSTRUMENTATION 

2.1  Flow  Apparatus 

3 

The  experiments  were  performed  in  the  Brown-Roshko 
flow  facility  designed  to  produce  a turbulent  shear  flow  between 
two  streams  of  different  gases. 

Basically  the  facility  consists  of  two  supply  lines,  each 
one  coming  from  eight  2000psi  bottles,  which  supply  two  gas 

i 

streams  that  are  brought  together  at  the  exit  of  two  nozzles 
which  have  been  contracted  down  to  1"  x 4"  each.  The  over- 
all size  of  the  test  section  is  2"  x 4"  x 12". 

The  settling  chambers  and  screens  in  the  supply  lines 
are  sufficient  to  reduce  the  free  stream  turbulence  level  in  the 
test  section  to  less  than  0.5%.  The  sidewalls  of  the  test  section 
are  adjustable  to  help  remove  any  pressure  gradients  in  the 
flow.  The  entire  test  chamber  can  be  pressurized  up  to  ten  at- 
mospheres in  order  to  vary  the  Reynolds  number.  Steady  state 
flow  is  usually  established  in  less  than  500  milliseconds  after 
the  flow  is  initiated. 

A traversing  mechanism,  which  incorporates  a digitally 
controlled  motor  to  step  the  probes  in  increments  of  .001", 
moves  the  probes  across  the  test  section  at  up  to  500  steps 
per  second. 

2.2  Instrumentation 

2.2a  Pressure  Data 

Two  types  of  pressure  data  were  required  for  this  in- 
vestigation: static  pressure  scans  along  the  test  section  used 
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in  minimizing  the  pressure  gradients  in  the  flow  direction 
(see  Section  3.1)  and  dynamic  pressure  profiles  across  the 
mixing  layer.  The  static  pressure  scans  were  accomplished 
with  the  use  of  a fast  pressure  scanner  (Scanivalve,  type- 
W1260),  a device  which  sequentially  communicates  each  of  the 
static  pressure  ports  along  both  sides  of  the  test  section  to 
the  input  of  a pressure  transducer  (a  Datametrics  electronic 
manometer,  1014A,  and  Barocel  differential  pressure  sensor, 
511-10).  The  reference  for  the  sensor  was  generally  the 
static  port  at  x = 1".  A pitot  tube,  connected  to  the  same 
type  of  pressure  transducer,  was  used  to  obtain  the  dynamic 
pressure  profiles. 

2.  2b  Concentration  Data 

In  order  to  determine  the  extent  of  mixing  within  the 
shear  layer  and  to  compute  the  mean  velocity  profiles  from  the 
dynamic  pressure  profiles,  the  local  composition  within  the 
layer  had  to  be  determined.  For  these  measurements  an  im- 
proved version  of  the  Brown-Rebollo  aspirating  probe^  was 
developed,  under  the  guidance  of  G.  L.  Brown,  to  meet  the  high 
frequency  and  small  spatial  resolution  requirements  (Fig.  2a). 

The  new  probe  is  estimated  to  have  a frequency  response 
of  100  kHz  and  a spatial  resolution  of  approximately  0.  1mm. 

The  signal  to  noise  ratio  was  sufficient  to  detect  concentration 
fluctuations  of  H.  A complete  account  of  the  probe  construction 
and  response  estimations  is  given  in  Appendix  A.  A standard 
constant  temperature  hot-wire  circuit  designed  and  built  at 
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Caltech  by  Anthony  Perry*  and  Brian  Cantwell**  was  used  to 
operate  the  hot  wire  in  the  probe. 

2.3  Data  Acquisition  System 

The  short  run  times  (the  Brown-Roshko  apparatus  is  a 
blow -down  type  facility  with  run  times  of  2-6  seconds)  and  the 
high  frequencies  of  the  experiments  necessitated  the  use  of  a 
high-speed  data  acquisition  system.  Therefore,  all  data  were 
taken  (and  reduced)  with  the  GALCIT  mobile  computer  con- 
trolled data  acquisition  system  (the  "Solo  System").  This  system 
has  the  capability  of  computer  controlled  data  acquisition,  re- 
duction, storage  and  output.  The  heart  of  the  system  is  a 
Hewlett-Packard  2100  computer  (32K  memory).  The  data 

acquisition  itself  is  accomplished  with  a sixteen  channel,  16  bit 

5 

analog-to-digital  converter  capable  of  5 x 10  single  channel 

5 

conversions  per  second  or  4 x 10  sequential  conversions  per 
second.  The  data,  before  or  after  reduction,  can  be  stored 
on  either  a Hewlett-Packard  read/write  magnetic  tape  deck,  7970E, 
or  a Hewlett-Packard  storage  disc,  7901  A.  The  output  devices 
with  the  system  are  a Teletype  and  a Hewlett-Packard  digital 
plotter,  72 10 A. 

The  required  interfacing  between  the  Solo  System  and  the 
experimental  apparatus  will  be  briefly  discussed  in  connection 
with  the  data  acquisition  procedures  in  the  next  section. 


* Reader,  University  of  Melbourne,  Australia 

**  Research  Fellow,  California  Institute  of  Technology 
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IIL  DATA  ACQUISITION  PROCEDURES 
3.1  Pressure  Gradient  Data 

One  of  the  boundary  conditions  for  the  experiments  was 
to  have  the  shear  flows  in  a zero  pressure  gradient  region. 
However,  any  vertical  fluid  column  has  a "natural"  hydrostatic 
pressure  gradient  equal  to  the  fluid  density  multiplied  by  the 
gravitational  acceleration.  The  effect  of  this  gradient  can  of 
course  be  neglected  in  experiments  in  which  the  fluid  is  air 
or  in  experiments  in  which  all  fluids  are  of  equal  densities  if 
care  is  taken  in  pressure  measurements. 

In  this  experiment,  which  had,  in  effect,  a column  of 
helium  in  direct  communication  with  a column  of  nitrogen, 
this  effect  could  not  be  neglected  or  eliminated,  only  minimized. 
(A  discussion  of  this  effect  and  the  minimization  procedure  are 
in  Appendix  B.)  In  order  to  minimize  this  effect,  the  static 
pressure  gradient  along  the  test  section  had  to  be  determined. 
This  was  done  by  sequentially  sampling  the  static  pressure 
ports  on  the  side  walls  of  the  test  section  with  the  use  of  a 
Scanivalve.  After  allowing  250  milliseconds  after  each  step 
of  the  Scanivalve  for  the  pressure  reading  to  reach  a steady 
state,  the  pressure  was  sampled  1000  times  in  250  milliseconds. 
The  readings  were  then  averaged  for  each  port  and  the  gradi- 
ent calculated.  The  ports  at  x = 2",  3",  4",  5",  6",  7"  and  8" 
on  one  wall  and  at  x = 4"  and  8"  on  the  other  were  sampled. 
The  side  walls  were  then  adjusted  to  give  the  optimal  gradient. 
The  side  wall  positions  were  different  for  each  free  stream 
velocity  setting. 


* 
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3.2  Traverse  Data 

The  stated  goal  of  determining  the  extent  of  mixing  in 
the  shear  layer  as  a function  of  cross-stream  position  and 
Reynolds  number  required  measuring  the  time  histories  of  the 
concentrations  of  the  gases  at  various  x and  y locations  with- 
in the  layer  for  a variety  of  free-stream  conditions.  The 
procedure  for  accomplishing  this  was  as  follows.  The  tank 
(ambient)  pressure  and  the  free-stream  velocities  were  ad- 
justed to  the  desired  settings.  The  concentration  probe  and 
a pitot  tube  both  mounted  approximately  one  inch  apart  on 
the  traversing  mechanism  were  placed  at  a selected  x station 
and  then  moved  across  the  shear  layer  alternately  sampling 
and  then  stepping. 

The  two  probes  were  sampled  sequentially  for  150 
milliseconds  at  each  point  at  a rate  of  80,000  readings  per 
second  for  a total  of  12,000  data  points.  They  were  then 
stepped  .075  inches  to  the  next  y location  (stepping  rate  was 
0.5  inches  per  second  for  150  milliseconds)  where  another 
12,000  data  points  were  taken.  After  each  data  acquisition 
period,  while  the  probes  were  moving,  the  data  were  trans- 
ferred from  computer  memory  to  disc,  thus  freeing  the  memory 
for  the  next  batch  of  data.  The  six  second  run  time  limitation 
allowed  a maximum  of  21  points  across  the  layer  (1.5  inches). 


After  each  six  second  run,  the  data  were  automatically  trans- 
ferred from  disc  to  magnetic  tape  and  then  back  into  computer 


memory  for  reduction.  The  reduced  data  were  then  stored  on 
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magnetic  tape  and  displayed  on  the  plotter.  Each  traverse 
was  then  repeated  four  times  and  the  mean  values  from  each 
traverse  were  averaged.  This  process  was  repeated  for  four 
x locations,  two  different  free-stream  velocity  settings  and  at 
ambient  pressures  of  one,  four  and  eight  atmospheres  for  a 
total  of  nearly  25  million  raw  data  points.  The  electronic 
interfacing  between  the  computer  and  the  experimental  appa- 
ratus consisted  of  a master  clock  (set  at  the  data  acquisition 
frequency,  usually  80  kHz  ) which  was  gated  alternately  to 
first  allow  the  12,000  pulses  to  go  to  the  A/D  converter  and 
then  divided  (to  500  Hz  ) and  gated  to  the  stepping  motor  for, 
usually,  75  pulses.  This  allowed  for  alternately  sampling  at 
80  kHz  and  then  stepping  at  500Hz. 

3.3  Power  Spectrum  Data 

In  order  to  reduce  the  time  histories  of  concentration 

to  power  spectral  density  curves  in  a reasonable  amount  of 

12  13 

time,  a fast  fourier  transform  routine  ' was  used.  This  re- 
quired 2n  data  points  taken  at  a frequency  at  least  twice  the 
maximum  frequency  of  interest  (always  less  than  25kHz  ex- 
cept for  one  run  at  Uj  = 3000cm/sec).  Therefore,  a data 
acquisition  rate  of  50kHz  was  used  with  a cut-off  filter, 
Krohn-Hite  3202A,  set  at  25kHz  on  the  output  of  the  concen- 
tration probe.  However,  the  computer  memory  limited  the 
number  of  data  points  which  could  be  reduced  to  4096  (~  80msec 
sampling  time).  This  is  not  sufficient  time,  however,  to  give  good 
resolution.  Therefore,  records  of  16,384  data  points  were  taken 
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IV.  DATA  REDUCTION 

The  first  step  in  reducing  the  data  acquired  by  the  preced- 
ing procedures  was  to  convert  the  voltage  output  of  the  concentra- 
tion probe  to  molar  concentrations.  This  procedure  is  explained 
in  Appendix  A. 

Discussed  in  the  following  sections  are  the  various  ways 
that  these  time  histories  of  concentration  and  dynamic  pressure 
were  converted  to  meaningful  results. 

4.  1 Mean  Profiles 

The  first  reduction  process  performed  on  the  data  was  to 

ensemble  average  the  data  for  each  point  across  the  mixing  layer 

Z~  2 — 

from  each  probe.  This  gave  PU  (y)/P[Uj  and  C(y).  The  con- 
centration profiles  were  then  converted  to  mean  density  profiles 
(p/pj  = C + (1.  - C)  p ^ /P  j )-  -An  estimation  of  the  velocity  pro- 
file was  found  by  dividing  the  mean  dynamic  pressure  by  the  mean 

density  and  taking  the  square  root.  This  procedure  neglects 

2 2 

several  fluctuation  correlations  (p'u'f  u'  , p'u'  ).  However, 

Rebollo^  showed  that  the  effect  is  negligible  (at  least  for  the 
2 2 

case  P2U£  /p^Uj  = 1). 

4.  2 Intermittency  and  Mean  Profiles  within  the  Turbulence 

Since  the  flow  had  a near  unity  Schmidt  number  (v/D),  the 
concentration  histories  could  be  used  to  estimate  the  intermittency 
function  across  the  mixing  layer  (percent  of  time  that  the  fluid  is 
turbulent).  If  the  concentration  at  some  x,  y,  and  t is  other  than 
0 or  1 , it  indicates  that  by  some  process  (ultimately  diffusion) 
the  fluids  from  the  two  sides  of  the  layer  have  become  molecularly 
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mixed.  However,  since  momentum  will  diffuse  approximately  the 
same  as  molecules  (Schmidt  number  unity)  this  also  indicates  that 
vorticity  will  have  diffused  into  the  sampled  region.  If  this  is  an 
indication  of  turbulence,  then  the  percent  of  time  that  the  concentra- 
tion is  other  than  1 or  0 will  give  an  estimation  of  the  intermit- 
tency  function  across  the  layer.  In  actual  fact  three  curves  can 
be  defined:  (1)  the  probability  of  finding  pure  gas  1 (within  3%), 

(2)  the  probability  of  finding  pure  gas  2,  and  (3)  intermittency 
(which  is  unity  minus  the  two  above  probabilities). 

This  procedure  of  defining  turbulent  and  nonturbulent  regions 
also  facilitated  estimating  mean  values  in  the  turbulent  regions 
(subscripted  T)  and  mean  values  in  the  nonturbulent  regions  (sub- 
scripted I).  Actually,  the  concentration  data  were  put  into  three 
different  classifications  as  illustrated  in  the  drawing.  Data  from 

(a)  were  used  to  find  means  in  nonturbulent  regions,  data  from 

(b)  were  used  to  find  mean  and  fluctuating  quantities  in  the  turbu- 
lent regions,  and  data  from  (a),  (b)  and  (c)  were  used  to  find 
overall  mean  and  fluctuating  quantities. 


Time 
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Thus,  the  following  quantities  were  measured  or  calculated: 

C(y)  - ensemble  average  of  concentration  histories  at 

each  y location. 

2 2 

PU  (y)/pjUj  - ensemble  average  of  dynamic  pressure 
5 histories  at  each  y location. 

U(y)/Uj  - mean  velocity  profiles. 

C(y)j  - average  of  concentration  data  in  the  nonturbulent 
regions  (data  within  3%  of  1 or  0). 

C(y)^,  - average  of  concentration  data  in  the  turbulent 
regions. 

Yj(y)  - probability  of  finding  pure  gas  1 (within  3%). 

Y2(y)  - probability  of  finding  pure  gas  2 
Y(y)  - Intermittency  ( 1 - - y2) 

( ),j,  - various  fluctuating  quantities  calculated  from  data 

in  the  turbulence  (as  explained  above).  These 
quantities  will  be  described  in  the  next  section. 

4.  3 "Unmixedness" 

One  of  the  goals  of  this  investigation,  to  determine  the 
extent  of  mixing  in  a turbulent  shear  layer,  required  a "yardstick" 
with  which  to  measure  the  extent  of  mixing.  The  definition  of  one 


such  "yardstick" 
where 


the  RMS  fluctuation  of  the  concentration, 


C'  = C - C . 


4.  1 


. 


However,  this  quantity  does  not  have  a well-defined  upper  bound 
and,  if  the  gases  in  the  layer  are  perfectly  unmixed  (concentra- 
tion data  are  always  either  1 or  0),  different  values  of  C 
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give  different  values  of  ^ . However,  for  comparison  with  other 
investigations,  this  quantity,  along  with  j,,  were  calculated 

where  C'  = C^,  - C^,,  C^,  being  concentration  data  within  the 
turbulence. 

A more  meaningful  measure  of  the  extent  of  mixing  comes 
from  chemical  reaction  rate  estimations.'  For  the  irreversible 
reaction 

A + B - P 4.  2 


where  A and  B are  reactants  and  P is  the  product,  the  following 

14 


conservation  equations  can  be  derived  : 

DC, 


Dt 


= - ki<t)Cacb 


4.  3 


DC, 


DT  = - K2<T>  CACB 


4.  4 


where  is  the  substantial  derivative,  IC  is  a reaction  rate,  C. 
is  the  concentration  of  species  i,  and  T is  temperature.  This 
approach  can  be  related  to  the  present  investigation  if  it  is  sup- 
posed that  A is  diluted  in  the  helium  and  B in  the  nitrogen. 

Then: 


CA  = CA  + CA  = CAo  |CHe>  * CAo  <CHe  + CHe» 


and 


4.  5 


CB  = CB  + C B = CBo  (1  - CH*>  ■ CB."  - CHe  * CH*>  • 


where  C.  is  the  free  stream  concentration  of  species  i. 
10 


* This  measure  of  mixedness  was  suggested  by  Coleman  duPont 
Donaldson. 
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DC 


Dt 


- Kl( 


A t r I r*  /> 


CA  CB  + CA  C'b- 


- K1  CA  C B 11  + 


c'ac'b 


4.  7 


A B 


- K1  CA  CB  11 


Cke2 


This  assumes  that  K'.  = 0.  Thus,  the  quantity 

— 2 — — 

1 - C'  /C  (1-C)  = M 


is  an  indication  of  how  well  the  fluids  are  mixed.  If  it  is  1(C'  = 0) 
then  the  fluids  are  perfectly  mixed  and  the  mean  values  (C^)  can  be 
used  to  estimate  reaction  rates.  On  the  other  hand,  if  this  quantity 
is  0,  the  fluids  are  completely  unmixed  (the  concentration  is  always 


0 or  1)  and  the  mean  values  give  a very  poor  estimation  of  reaction 

~2 


rates.  Therefore,  the  quantity  C'  /(1-C)  C = 1-M  (the  unmixedness 
of  the  fluids)  was  calculated. 


17 


This  definition  of  unmixedness  will  obviously  amplify  the 
experimental  "noise"  in  measuring  C'  as  C -*  1 or  0.  Therefore, 
another  measure  of  unmixedness  which  avoids  this  problem  and 
still  has  the  proper  limits  was  defined. 

This  measure  of  unmixedness  is  defined  as 


where 


J(C- C)dtj  + J (C-C)dt2 


UM  = 


J11*'  d‘>  + J 


(C)dt, 


tj  = time  when  C > C 
t 2 = time  when  C < C 


4.  8 


This  integral  definition  is  illustrated  in  the  following  drawing. 


UM  is  actually  the  sum  of  the  crosshatched  areas  divided  by  the 
sum  of  the  areas  of  the  rectangles.  Obviously,  if  the  fluids  are 
perfectly  mixed  (C(t)  = U)  UM  = 0,  and  if  the  fluids  are  perfectly 
unmixed  (C(t)  = 1 or  0),  UM  = 1.  Therefore,  UM  has  the  proper 
limits.  Furthermore,  it  does  not  present  any  experimental 
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problems  for  C -*  1 or  0.  As  will  be  seen  later  in  the  discussion 
of  a chemically  reacting  shear  layer,  this  definition  also  has  an 
interpretation  relating  the  amount  of  product  formed  to  the  amounts 
of  reactants  entrained. 

4.  4 Determination  of  Dividing  Streamline 

In  order  to  properly  show  all  of  the  traverses  for  the  var- 
ious x locations  in  nondimen sional  form,  a common  reference  point 
for  all  of  the  profiles  has  to  be  determined.  This  point  is  com- 
monly taken  to  be  the  intersection  of  the  dividing  streamline  with  the 
y-axis  of  the  data  traverse.  This  point  is  found  by  combining  the  con- 
tinuity and  momentum  equations  to  derive  the  following  equation: 

y* 


d 


/ ^ (u,  ' ’)  dy  = £ J pP,\  («i  ' “,) 


dy 


4.  9 


where  ( )j  denotes  free  stream  conditions  on  side  1,  ( denotes 
free  stream  conditions  on  side  2 and  y*  is  the  location  of  the 
dividing  streamline  intersection  point.  The  application  of  Eq.  4.  9 
to  find  y*  for  a shear  layer  having  similarity  properties  is  rather 
straightforward.  However,  in  applying  it  to  an  inhomogeneous 
wake  (U£  = Uj,  p^  t P j ) more  care  must  be  taken.  One  possible 
approach  is  as  follows.  Let 
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■ ~ I ('  ■ ^)dv'  *2= kR) 


dy 


4.  10 


Then,  combining  Eqs.  4.9  and  4.  10  gives 


Ideas  similar  to  these  were  suggested  by  Dr.  Erik  Storm  of 
Lawrence  Livermore  Laboratory. 
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a.  , . /M*  jL  e 

dx  W1  \U2f  dx  °2  * 

£ 01  = - 4 h-  <U2  - Ul» 


8j  can  be  interpreted  as  the  momentum  delect  on  side  1 and  8^ 
the  momentum  defect  on  side  2.  Equation  4.  12  states  that  if 
9j  increases  (decreases),  0£  must  decrease  (increase)  by  exactly 
the  same  amount.  In  other  words,  6^  + 0^  = 9 = constant,  as 

must  be  the  case  (the  momentum  defect  for  a wake  is  constant). 

If  8 j and  9^  are  continuously  changing,  then  one  side  of  the  wake 
must  be  continuously  gaining  momentum  and  the  other  side  continu- 
ously losing  momentum.  As  this  situation  can  not  continue  indef- 


initely, Eq.  4.  12  must  be  interpreted  as 


dx  2 


0 


or  0 j /8  2 = const.  4.13 

Since  an  estimate  of  the  constant  in  Eq.  4.  13  can  be  made  from 
the  initial  boundary  layers  on  the  splitter  plate,  Eq.  4.13  together 
with  Eq.  4.  10  can  be  used  to  find  y*  for  the  inhomogeneous  wake. 
(This  estimation  of  y*  neglects  all  fluctuation  correlations.  ) 

4.  5 Entrainment 

With  today's  large  structure  view  of  turbulence,  entrainment 
becomes  a difficult  quantity  to  define.  The  traditional  definition 
of  entrainment  for  a shear  layer. 


. 

r 


trarily  chosen  limit,  raises  two  questions,  1.  How  should  y. 


logically  be  chosen?  2.  Should  the  fluid  which  is  within  the  bound  of 
yj  and  y^,  but  has  not  yet  become  part  of  the  turbulence,  be  part  of  the 
entrainment  calculation?  If  the  entrainment  rate  is  to  provide 
information  as  to  the  proportions  and  rate  in  which  the  two  fluids 
on  either  side  are  entering  the  mixing  layer  and  becoming 
"entangled"  or  mixed,  then  a suitable  definition  is  needed.  One 
such  definition  is 


where  Y(y)  is  the  intermittency  function  and  (PU)^,  is  the  average 

within  the  turbulence  of  the  product  of  density  and  velocity.  The 
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complete  derivation  of  this  definition  is  given  by  Brown  . For 
this  investigation  (pU)j  was  approximated  by  p,p  U. 

It  is  not  entirely  clear  whether  Eq.  4.  14  or  Eq.  4.  15  should 
be  used  to  describe  the  entrainment.  Both  definitions  describe  the 


amount  of  fluid  which  has  crossed  an  arbitrary  boundary.  The 
limit  used  in  Equation  4.  14  is  an  overall  boundary  of  the  shear 
layer  based  on  mean  quantities  and  the  limit  used  in  Equation  4.  1 5 
is  the  boundary  between  the  turbulent  and  nonturbulent  regions. 
However,  for  clarity,  the  quantities  calculated  using  Equation  4.15 
will  be  referred  to  simply  as  the  entrainment  rates  throughout  the 
rest  of  this  work.  This  is  done  since  these  rates  are  more  mean- 
ingful in  relation  to  the  chemical  analysis  presented  in  Section  7 
than  the  rates  calculated  using  Equation  4.  14.  Further,  entrain- 
ment into  the  turbulence  is  more  in  accord  with  the  large  structure 
view  of  turbulence  than  is  simply  a gross  measure  of  entrainment. 

The  normalization  was  chosen  so  that  entrainment  rates  can  be  com- 
pared directly  without  having  to  take  differences  in  density  into 
account. 

4.  6 Power  Spectrum 

As  discussed  in  Section  3.  3,  four  records  of  4096  concentra- 
tion readings  each  were  taken  at  a rate  of  50  kHz  for  various  x and  y 
locations  and  various  free-stream  velocities.  These  four  records, 
after  being  put  through  a Hanning  window,  were  reduced  to  four 
spectra  using  a fast  fourier  transform  routine.  These  were  then 
averaged  to  produce  a single  spectrum.  This  curve  was  then  further 
reduced  by  averaging  adjacent  points  within  the  curve.  The  averaging 
window  (number  of  adjacent  points  to  be  averaged)  was  exponentially 
increased  from  one  for  the  low  frequency  points  to  approximately  30 
for  the  highest  frequencies. 
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V.  RESULTS  AND  DISCUSSION 

The  initial  part  of  this  investigation  was  to  measure  the 
extent  of  mixing  for  a shear  layer  of  fixed  parameters  for  as 
wide  a Reynolds  number  range  as  possible.  Following  are  the 
parameters  which  were  kept  constant: 

(a)  A two-dimensional  shear  layer  between  high- 

P? 

speed  helium  and  low-speed  nitrogen  (s  = — — = 7.0) 

P1 

(b)  p(x)  = constant  (p  being  static  pressure) 
and  (c)  r = U2/U1  = 0.38  . 

Following  are  the  parameters  which  were  changed  in  order  to 
vary  the  Reynolds  number: 

(a)  Uj  (and  thus  U2):  1000,  2000  and,  for  Pq  = 1 atm., 
3000cm/  sec. 

(b)  Pq  (pj  and  P2)  : 1,  4 and  8 atmospheres 
and  (c)  x:  2,  3,  4 and  5 inches. 

After  these  data  were  taken,  two  other  cases  were  investi- 
gated. In  order  to  determine  what,  if  any,  effect  the  inhomogeneity 
of  the  flow  has  on  the  mixing,  the  mixing  in  a constant  density 
shear  layer  was  measured  (r  = 0.38,  s = 1.0)  and  compared  to 
the  previous  results.  For  this  case,  the  high-speed  gas  was  ni- 
trogen and  the  low -speed  gas  a mixture  of  argon  and  helium 
arranged  to  have  the  same  density  as  nitrogen.  (See  Appendix 
C for  a discussion  of  the  method  used  to  measure  the  density 
of  the  gas  mixture.)  The  concentration  probe  was  able  to  re- 
solve the  concentration  of  nitrogen  in  this  mixture  to  within  4%. 

For  this  case  Pq  = 4 atms,  Uj  = 1000  cm/sec,  and  x = 4 inches. 
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The  third  case  investigated  was  a wake,  but  not  a wake 
in  the  usual  sense  of  the  word.  This  case  was  actually  the 
investigation  of  the  wake  produced  by  the  splitter  plate  which 
initially  separates  the  two  fluids.  The  flow  was  obtained  simply 
by  settling  U2  equal  to  Uj  (r  = 1.0,  s = 7.0).  This  case  was 
investigated  to  compare  the  mixing  in  a wake  to  that  of  a shear 
layer.  Since  the  structure  of  a wake  is  quite  different  from 
that  of  a shear  layer,  the  mixing  process  may  be  different. 

For  this  case  Uj  = = 2000cm/sec.,  Pq  = 1 atm.,  and  x = 

2,  3,  4,  5 and  6 inches. 

In  the  following  sections  are  the  results  for,  first  of  all, 
the  inhomogeneous  shear  layer,  then  the  homogeneous  shear 
layer,  and,  finally,  the  wake. 

5.1  Shear  Layer  (r  = 0.38,  s = 7.0) 

5.  la  Mean  Profiles  and  Spreading  Rates 

First  of  all,  in  order  to  acquire  a "feel"  for  the  extent 
of  the  fluctuations  in  the  concentration  of  helium  within  the 
shear  layer,  the  probe  was  continuously  traversed  across  the 
layer  from  the  nitrogen  to  the  helium.  The  result  is  shown  in  Fig- 
ure 3a.  Similar  results  were  obtained  by  Brown  and  Roshko. 
Several  questions  were  immediately  raised  by  this  curve.  Why 
do  the  fluctuations  appear  to  be  larger  on  the  nitrogen  side  than 
on  the  helium  side?  Why  does  the  concentration  appear  to 
spend  a large  amount  of  the  time  near  C = 0.8?  The  answers 
to  these  questions  are  discussed  in  the  section  on  entrainment 
rates.  As  observed  by  Brown  and  Roshko,  shadowgraphs  of  the 


flow  at  various  Reynolds  numbers  (Plate  la,  c,  e,  g ) show  that 
large  structures  exist  in  the  shear  layer  and  these  structures 
appear  to  be  entraining  the  pure  free-stream  fluid  in  such  a 
way  that  it  can  exist  in  an  unmixed  state  across  the  entire 
layer.  Thus,  large  fluctuations  in  concentration  are  expected. 
This  will  be  further  discussed  in  Section  6. 

Figure  3b  shows  the  mean  concentration  profiles  as  de- 
scribed in  Sections  3.2  and  4.2.  A characteristic  length  can 
be  defined,  and  thus  a growth  rate  can  be  found,  from  this  curve. 
If  the  "edges"  of  the  layer  are  defined  as  the  points  C = 0.1 


and  C = 0.9  then 


6 = 0.147  (x  - x ) 

c o 


Brown  and  Roshko  used  C = 0.01  and  C = 0.99  as  the  edges  of 
the  shear  layer  and,  thus,  their  measured  growth  rate  was  some- 
what larger. 

Figures  3c  and  3d  are  the  mean  dynamic  pressure  and 
the  mean  velocity  profiles,  respectively.  Another  length,  one 
which  can  also  be  defined  for  homogeneous  as  well  as  inhomo- 
geneous shear  layers,  can  be  found  from  Figure  3d.  This 
length  is  the  maximum  slope  thickness  of  the  velocity  profile 
or  the  vorticity  thickness  and  is  measured  to  be 

6uu  = °-119  (x  - xq)  5.2 

This  characteristic  length  will  be  used  throughout  the  rest  of 
this  work  and  will  be  labeled  6. 

Shown  in  Figure  4 are  a few  typical  time  histories  of  the 


concentration  taken  at  various  T|  locations  in  the  layer.  L is 
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the  approximate  size  of  the  large  structure  for  that  x location. 

It  is  important  to  note  that  these  mean  profiles  were  un- 
affected by  varying  the  Reynolds  number. 

5.  lb  RMS  Concentration  Fluctuations 

A rather  surprising  result  was  found  when  the  fluctuating 

quantities  were  compared  on  similarity  plots.  Two  curves  were 

found  for  each  quantity  plotted,  one  for  the  profiles  taken  at 

Reynolds  numbers  less  than  2x10  (Re  = &AU/V  ) and  one  for 

1N2 

c 4 

the  profiles  taken  at  Reynolds  numbers  greater  than  2x10 

This  Reynolds  number  will  henceforth  be  called  the  "critical" 

Reynolds  number  (Rec).  Data  curves  above  Rec  will  be  labeled 

Re  + and  those  below  will  be  labeled  Re  . This  effect  is 
c / , . . c 

shown  for  { C 2 and  (y  C 2 ) T in  Figures  5a  and  b.  The  Re^ 
data  are  shown  by  the  dashed  curves.  The  fluctuations  in  con- 
centration are  always  less  above  the  critical  Reynolds  number 

than  below.  Data  taken  near  Re  fell  between  the  two  curves. 

c 

One  thing  to  be  noted  in  Figure  5a  is  that  the  nitrogen  side 
appears  to  have  a higher  degree  of  fluctuations  than  the  helium 
side.  This  will  be  partially  explained  in  the  discussion  on  en- 
trainment. Two  things  should  be  noted  about  Figure  5b:  the 
profile  is  flatter  than  that  of  Figure  5a  and  the  Reynolds  number 
effect  is  more  pronounced  for  the  turbulent  zone  calculation 
than  for  the  calculation  for  the  entire  data  record  (hereafter 
called  the  normal  calculation).  These  observations  will  be  use- 
ful in  determining  why  there  is  a Reynolds  number  effect  in 
the  first  place.  This  will  be  discussed  in  Section  6. 
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5.  lc  Intermittency  and  Entrainment 

This  same  Reynolds  number  effect  can  be  observed  in  the 
intermittency  similarity  curves.  Figure  5c- 1 is  the  probability  of  find- 
ing pure  nitrogen,  Figure  5c-2  is  the  probability  of  finding  pure  helium 
and  Figure  5c-3  is  the  intermittency  function.  The  dashed  curves  are 
for  the  Rec  data.  For  the  Rec+  data,  the  scatter  may  be  judged  from 
the  plotted  points;  for  all  the  Rec~  data  the  scatter  was  smaller  by 
about  one  half.  It  can  be  seen  from  Figure  5c- 1 and  2 that  below  Rec> 
the  probability  of  finding  pure  gas  is  higher  across  the  entire  shear 
layer.  Thus  the  intermittency  function  for  Rec  is  more  narrow  and 
reaches  a lower  maximum  value  than  for  Rec+.  Therefore,  the  pure 
fluid  which  is  being  engulfed  into  the  shear  layer  is  being  entrained 
into  the  turbulence  (being  mixed  molecularly)  faster  above  Rec  than 
below.  Pure  fluid  extends  across  most  of  the  shear  layer  less  often 
above  Rec  than  below.  Besides  the  intermittency  function,  the  entrain- 
ment rates  into  the  turbulence  are  other  quantitative  measures  of 
this  effect.  Using  Equation  4.  15  the  entrainment  rates  were  found 

to  be  Ve.  /U,  = 0.  062  and  Ve,/U,  = 0.019  for  Re  " and  Ve.  /U.  = 

it  cl  c 11 

0.067  and  Ve^Uj  = 0.020  for  Rec+.  This  is  approximately  an 
eight  percent  increase  in  the  entrainment  rates  into  the  turbulence. 
The  traditional  measurement  of  entrainment  (Eq.  4.  14)  would  not 
be  affected  by  varying  the  Reynolds  number  as  the  similarity  pro- 
files of  the  density  and  velocity  are  the  same  above  and  below  Rec. 
The  significance  of  the  Reynolds  number  affecting  the  fluctuating 
quantities  but  not  the  mean  quantities  will  be  discussed  in  Section  6. 
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Besides  these  increases  in  the  entrainment  rates,  another 
important  point  to  note  is  that  over  three  times  as  much  helium 
is  being  entrained  as  nitrogen.  This  fact  accounts  for  several 
of  the  points  brought  out  in  earlier  discussion.  Due  to  this 
difference  in  entrainment  rates,  the  large  structures  will  have 
a mean  concentration  near  0.8.  Thus  the  output  from  the  con- 
centration probe  will,  for  a large  percentage  of  the  time,  be 
near  this  value  (as  seen  in  Figs.  3a,  3b  and  4).  This  also 
partially  explains  why  the  fluctuations  are  larger  on  the  nitrogen 
side.  On  that  side,  as  the  large  structures  pass  by  the  probe 
(see  drawing)  the  concentration  will  fluctuate  from  approximately 
0.8  to  0.0.  On  the  helium  side,  however,  the  concentration  will 
fluctuate  from  0.8  to  1.0.  This  only  partially  explains  the  larger 
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fluctuations  on  the  nitrogen  side  because,  as  can  be  seen  from 
Fig.  5b,  the  turbulent  zone  calculations  still  give  larger  fluctu- 
ations on  the  nitrogen  side  than  on  the  helium  side,  although 
the  difference  is  not  as  great  as  in  the  normal  calculations 
(Fig.  5a). 

5.  Id  Unmixedness 

Figure  6 shows  the  results  of  the  unmixedness  calculations. 
Figure  6a  shows  the  results  for  the  integral  definition  (Eq.  4.8) 
and  Figure  6c  shows  the  results  for  the  chemical  reaction  rate  defi- 
nition. Figures  6b  and  6d  are  the  turbulent  zone  calculations  for 
the  unmixedness.  As  with  the  RMS  calculations,  the  unmixedness 
values  are  always  less  above  Rec  than  below.  Also,  lower  values 
are  obtained  and  the  Reynolds  number  effect  is  more  pronounced 
for  the  turbulent  zone  calculations  than  for  the  normal  calculations. 

Several  qualitative  statements  can  be  made  from  the  data 
shown  on  these  curves.  Firstly,  from  Figure  6a,  fluid  which  has 
been  entrained  in  the  usual  sense  (Eq.  4.14)  is  still  in  a rela- 
tively high  state  of  unmixedness,  has  not  become  mixed  on  the 
molecular  level,  especially  on  the  nitrogen  side  and  below  Rec  . 

If  the  nitrogen  and  helium  were  reactive  with  each  other,  much 
of  the  fluid  within  the  layer  would  not  react  at  the  point  at 
which  it  has  been  entrained;  it  would  "swirl"  around  within  the 
layer  without  having  the  opportunity  to  mix  on  the  molecular 
level,  and  thus  react.  Of  course,  all  of  the  fluid  which  is  en- 
trained will  eventually  react,  but  only  a small  portion  of  the 
total  at  any  r location  will  have  been  reacted.  Even  the  fluid 


which  has  been  entrained  into  the  turbulence  is  in  a high  state 
of  unmixedness  (Fig.  6b). 

Since  information  about  the  degree  of  unmixedness  is  not 
usually  known. 
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is  sometimes  used  to  calculate  reaction  rates  instead  of  Eq.  4.  7 
It  can  be  seen  from  Figure  6c  that  this  calculation  would  be  in 
error  by  at  least  20%. 

These  qualitative  results  will  be  made  more  quantitative 
in  Section  7 where  these  results  are  used  to  predict  amounts 
of  product  in  a chemically  reacting  shear  layer. 

5.  le  Probability  Density  Function 

A useful  way  of  presenting  the  data  in  the  time  histories 
of  the  concentration  is  with  probability  density  functions  (P.D.F.'  s). 
P.D.F.'s  are  simply  graphs  of  the  percentage  of  time  that  the 
concentration  spends  between  two  limits  (normalized  to  give 
unity  area  under  the  curve).  Shown  in  Figure  7 are  the  P.D.F.'s 
for  various  p locations  across  the  layer.  All  values  above  8.0 
are  not  shown.  It  is  rather  interesting  to  note  that  all  of  the 
results  concerning  concentration  presented  thus  far  can  be  cal- 
culated from  the  P.D.F.'s.  That  is  one  of  the  reasons  why  at 
present  a considerable  amount  of  effort  is  being  spent  on  trying 
to  model  the  P.D.F.'s  of  scalar  quantities  in  turbulent  shear 
flows.  One  popular  practice  is  to  model  them  with  a Gaussian 
distribution.  As  can  be  seen,  at  no  point  in  this  case  do  the 
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P.D.F.'s  approach  a Gaussian  distribution. 

The  effect  of  the  difference  in  entrainment  rates  can  again 
easily  be  seen.  The  probability  of  finding  the  concentration  be- 
tween 0.8  and  0.9  is  relatively  high.  The  probability  of  finding 
pure  gas  across  the  layer  is  shown  by  the  spikes  on  the  ends 
of  the  curves.  This  particular  set  of  curves  is  the  average  of 
four  runs  at  x = 4 inches,  Uj  = lOOOcm/sec  and  Pq  = 4 atms. 

No  attempt  was  made  to  average  all  of  the  P.D.F.'s  above  Rec 
and  all  below  Rec  because  each  set  of  runs  measured  the  P.D.F.'s 
for  different  R locations.  These  P.D.F.'s  will  be  used  exten- 
sively in  Section  7. 

5.  If  Power  Spectral  Density  Curves 

In  order  to  gain  more  insight  into  the  observed  Reynolds 

number  effect,  spectra  of  the  concentration  fluctuations  were 

obtained  both  above  and  below  Re  for  various  p locations. 

c 

Shown  in  Figure  8a  are  curves  obtained  at  P = 0.0.  Curves 

marked  (1)  are  for  data  below  Re  , curves  marked  (2)  are 

c 

for  data  near  Rec  and  curves  marked  (3)  are  for  data  above 
Rec  . Curves  (1)  and  (3)  show  a marked  difference  in  high 
frequency,  small  scale  mixing  that  takes  place  above  and  below 
Rec;  at  f6/AU  = 10  the  two  sets  of  spectral  curves  are  separated 
by  more  than  two  orders  of  magnitude.  The  spectrum  above  the 
critical  Reynolds  number  is  qualitatively  similar  to  most  high 
Reynolds  number  turbulent  velocity  spectra.  It  is  significant 
that  the  low  frequency  end  of  the  spectrum  is  the  same  above 
and  below  the  critical  Reynolds  number.  Thus,  the  overall 
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mechanics  and  development  of  the  large  structure  is  essentially 
unaltered  by  the  increase  in  fine  scale  mixing  above  Rec  . 

The  differences  in  the  spectrum  for  two  "H  locations 
at  the  same  Reynolds  number  is  shown  in  Figure  8b.  This  com- 
parison is  indicative  of  the  difference  in  RMS  fluctuations  for 
the  two  locations. 

5.2  Shear  Layer  (r  = 0.38,  s = 1.0) 

5. 2a  Mean  Profiles.  Intermittency  and  Entrainment 
As  with  the  nitrogen-helium  shear  layer  in  order  to  get 
a feel  for  the  extent  of  unmixedness  of  this  constant  density 
shear  layer,  the  concentration  probe  was  continuously  sampled 
while  being  traversed  across  the  layer.  The  results  are  shown 
in  Figure  9a.  Again,  the  extent  of  the  fluctuations  of  the  concen- 
tration of  nitrogen  is  quite  large,  but,  in  this  case,  the  fluctuations 
are  rather  symmetric  about  the  middle  of  the  layer.  This  symmetry 
can  more  easily  be  seen  in  the  mean  concentration  profiles  shown 
i in  Figure  9b.  Shadowgraphs  of  this  flow  are  shown  in  Plate  2. 

4 

The  remarkably  flat  profile  of  the  turbulent  zone  mean  concen- 
tration is  noteworthy.  Similar  results  were  obtained  by  Fiedler1^ 
for  the  temperature  profiles  in  a slightly  heated  shear  layer. 

The  intermittency  and  mean  velocity  profiles  are  shown  in  Figures 
9c  and  d,  respectively.  The  entrainment  rates  calculated  for  this 
case  are  VCj/Uj  = 0.056  and  Ve2/Uj  = 0.042.  Thus,  only 
slightly  more  nitrogen  was  being  entrained  than  argon-helium 
mixture.  A few  of  the  time  histories  used  in  these  calculations 
are  shown  in  Figure  10. 
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5.  2b  Unmixedness  and  RMS  Fluctuations 

Due  to  instrumentation  problems  only  one  set  of  data  was  taken 
4 

for  this  case  (Re  = 2x10  ).  As  will  be  shown  in  Section  6,  this  is  above 
the  corresponding  critical  Reynolds  number.  The  unmixedness  profiles 
for  this  constant  density  shear  layer  are  shown  in  Figures  11a  and  b. 
Like  the  nitrogen  helium  shear  layer,  these  values  are  relatively  large 
and  always  less  for  the  turbulent  zone  calculations  than  for  the  normal 
calculations.  The  similarities  in  the  profiles  for  the  two  cases  (a  min- 
imum in  the  middle  with  larger  maxima  on  the  edges)  along  with  the 
fact  that  the  observed  large  structures  for  both  cases  are  quite  similar 
may  indicate  that  the  mixing  processes  involved  are  quite  similar.  In 
other  words,  the  mixing  process  may  be  much  more  a function  of 
r (t^/Uj)  than  s (P^/P^).  The  RMS  profiles  are  shown  in  Figure  11c. 

5.  2c  Probability  Density  Function 


The  P.  D.  F.  's  for  this  case  are  shown  in  Figure  12.  The 
slight  difference  in  entrainment  rates  is  evidenced  by  the  slight  skew- 
ness of  the  P.  D.  F.  's  toward  the  nitrogen  side.  They  are  otherwise 
apparently  symmetric  both  in  r|  and  C.  Again,  a Gaussian  fit  to  these 
P.  D.  F.  's  would  be  hard  to  justify. 

5.  3 Wake  (r  = 1.  0.  s = 7.  0) 


5.  3a  Mean  Profiles  and  Spreading  Rates 


The  continuous  traverse  profile  for  this  inhomogeneous 
wake  is  shown  in  Figure  13a.  Even  though  the  0 to  1 fluctuations  are 
not  as  pronounced  as  in  the  previous  two  cases,  the  concentration 
fluctuations  are  still  large.  The  profile  also  appears  to  be  surprisingly 
symmetric. 


The  mean  concentration  and  velocity  profiles  are  shown  in 

Figures  13b  and  c.  U is  the  maximum  velocity  defect.  The 

s 

scatter  in  the  velocity  profile  is  mainly  due  to  the  normaliza- 
tion. The  scatter  represents  less  than  1%  of  the  free-stream 
velocity,  but  about  10%  of  Us  . This  is  only  a rough  estimation 
of  the  actual  velocity  profile  since  this  calculation  neglects  all 
fluctuation  correlations.  For  this  case,  several  of  these  cor- 
relations (especially  P'u')  probably  can  not  be  neglected.  However, 
no  estimation  of  their  importance  can  be  made  from  these  data. 

A shadowgraph  of  a wake  at  this  Reynolds  number  is  shown  in 
Plate  3e. 

Again,  several  lengths  can  be  defined  from  these  curves. 

If  the  mean  concentration  curve  is  used  to  define  the  growth 
rate  and  the  "edges"  of  the  wake  are  defined  as  C = 0.1  and 
C = 0.9  then  for  0=9^  + 0^  = 0.  1 1 inches, 

,0.5 

(6c/0)  = 0.67((x  - xo)/0)  . 5.3 

If  the  mean  velocity  profile  is  used  to  define  the  growth  rate  of 
the  wake  and  the  "edges"  are  defined  as  the  intersections  of  the 
maximum  slope  lines  of  the  velocity  profile  with  the  line 

(U,  - U)/U  = 0 then 

8 / \°. 5 

(fiy/0)  = 1.06((x  - xo)/0)  5.4 

The  0.5  exponent  was  suggested  by  similarity  arguments  for  a 
homogeneous  wake  and  these  data  closely  fit  this  growth  law. 

Typical  time  histories  of  the  concentration  for  various  tj 
locations  are  shown  in  Figure  14. 
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5.  3b  RMS  Fluctuations,  Intermittency  and  Entrainment 
The  RMS  fluctuation  curves  are  shown  in  Figures  15a  and 
b.  In  contrast  to  the  results  for  the  previous  cases,  these  profiles 
are  relatively  flat.  The  remarkably  flat  profile  indicates 

that  within  the  turbulence  the  level  of  fluctuation  is  rather  constant 
across  the  wake.  The  absolute  level  of  fluctuations  is  slightly 
less  than  that  obtained  for  the  other  two  cases.  In  contrast  to  the 
case  for  the  shear  layer,  the  intermittency  (Fig.  15c- 3)  is  rather 
skewed  to  the  helium  side.  This  is  probably  due  to  the  difference 
in  entrainment  rates  for  the  two  fluids.  These  entrainment  rates 
were  calculated  to  be 


0.  5 


Vej/Uj 


and 


; °- 22 

v.-j/u,  = 0. 16  (xhr) 


5.  5 


Thus,  for  small  x,  the  wake  entrains  fluid  at  a faster  rate  than 
the  shear  layer. 

Other  than  this  skewness  to  the  helium  side,  this  inter- 
mittency function  is  quite  similar  to  those  obtained  for  the  other 
two  cases.  This  is  rather  remarkable  considering  the  difference 
in  the  large  structure  for  the  wake  and  the  shear  layer. 

5.  3c  Unmixedness 

The  unmixedness  profiles  for  the  wake  are  shown  in  Figure 
16.  Again,  the  profiles  are  remarkably  flat  and  the  values  of 
unmixedness  are  lower  than  for  the  shear  layer  (the  fluids  are 
more  mixed  on  a molecular  level).  Since  the  unmixedness 


& 


35 

profiles  for  the  wake  are  dissimilar  to  the  profiles  for  the  shear 
layer,  the  mixing  process  involved  is  probably  different  for  the  two 
cases.  This  of  course  is  expected  considering  the  difference  in 
the  large  structures  for  the  shear  layer  and  the  wake.  One  con- 
clusion which  can  be  drawn  from  these  unmixedness  profiles  and 
the  entrainment  rates  is  that,  for  small  x,  more  product  will  be 
produced  in  the  wake  than  in  the  shear  layer  if  the  two  fluids  are 
reacting.  This  is  because  more  molecular  mixing  and  more 
entrainment  are  occurring  in  the  wake  than  in  the  shear  layer. 

Of  course,  for  large  x,  the  entrainment  for  the  shear  layer  is 
much  larger  than  for  the  wake  and  thus  more  product  is  being  pro- 
duced in  the  shear  layer.  What  is  meant  by  "large"  and  "small"  x 
will  be  shown  in  Section  7. 

5.  3d  Probability  Density  Function 

Typical  P.  D.  F. 's  for  the  inhomogeneous  wake  are  shown  in 
Figure  17.  The  P.  D.  F. 's  are  quite  symmetric  both  in  T|  and  C. 
Again,  as  with  the  previous  P.  D.  F's,  a Gaussian  fit  would  be  hard 
to  justify.  As  mentioned  earlier,  these  results  will  be  made  more 
quantitative  in  Section  7 by  applying  the  data  to  a case  in  which 
the  amount  of  product  produced  in  a reacting  flow  can  be  predicted 
given  the  behavior  of  scalar  contaminants  in  a similar  nonreacting 
flow. 

Discussed  in  the  next  section  is  a photographic  study  under- 
taken to  gain  more  understanding  of  the  Reynolds  number  effect 
observed  in  the  inhomogeneous  shear  layer.  The  other  two  cases 
were  also  studied  to  see  if  any  evidence  for  a similar  effect  in 
each  case  could  be  found. 
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VI.  PHOTOGRAPHIC  INVESTIGATION  AND  REYNOLDS  NUMBER 
EFFECT 

4 

A close  look  at  the  shadowgraphs  of  the  nitrogen-helium 
shear  layer  at  various  Reynolds  numbers  reveals  why  this  observed 
Reynolds  number  effect  was  somewhat  expected.  Shown  in  Plate 
la,  c,  e and  g are  shadowgraphs  of  this  flow  at  various  Reynolds 
numbers.  The  critical  Reynolds  number  observed  in  the  measure- 
ments occurs  at  the  streamwise  position  approximately  one-third 
of  the  total  length  along  Plates  lc  and  le.  (Henceforth,  stream- 
wise  positions  will  be  labeled  as  the  fraction  of  the  total  length 
of  the  plate  in  question,  i.  e,  , this  position  will  be  labeled  x = 1/3.) 
It  is  quite  easy  to  discern  a marked  change  in  the  qualitative 
appearance  of  the  flow  at  this  point.  The  large  structures 
upstream  of  this  point  and  at  Reynolds  numbers  below  Rec  are 
still  evidently  quite  sharply  defined,  while  downstream  of  this  \ 

point  and  at  Reynolds  numbers  above  Re^  the  structures  exhibit 
much  more  three-dimensional  fine  scale.  The  sharpness  of  this 
transition  suggests  that  this  observed  Reynolds  number  effect  is 
caused  by  an  instability  of  some  kind.  Since  the  overall  behavior 
and  development  of  the  large  structures  are  not  greatly  affected 
by  this  instability  (the  linear  growth  rate  continuous  through  this 

transition  and  the  observation  of  large  structures  persists  at 

17 

Reynolds  numbers  much  above  Re^.  ),  it  was  felt  that  the  insta- 
bility was  associated  with  the  large  structure  itself  rather  than 
with  the  entire  flow  field.  In  order  to  reinforce  these  ideas,  a 
photographic  study  of  both  the  plan  and  side  views  of  the  shear 
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layer  was  undertaken.  The  results  of  this  study  are  presented  in 
the  next  section  and  discussed  with  regard  to  the  Reynolds  number 
effect  in  Sections  6.  2 and  6.  3. 

6.  1 Photographic  Investigation 

The  simultaneous  plan  and  side  views  for  this  investigation 
were  obtained  by  using  a beam  splitter  in  the  incident  parallel 
light  beam.  Only  half  of  the  plan  view  width  is  shown. 

Possibly  the  best  way  to  understand  what  is  happening  as 
the  shear  layer  develops  and  thus  to  see  what  is  causing  the 
Reynolds  number  effect  is  to  follow  the  development  of  the  shear 
layer  photographically  as  the  Reynolds  number  is  increased  past 
Rer.  (A  good  Reynolds  number  survey  of  the  shear  layer  is 
shown  in  Fig.  20  of  Brown  and  Roshko.  ) Shown  in  Plate  la  and  b 
is  the  shear  layer  at  a relatively  low  Reynolds  number.  The 
completely  two-dimensional  nature  of  these  large  structures  is 
shown  by  the  upstream  portion  of  both  the  plan  and  side  views. 
This  initial  roll-up  is  essentially  due  to  a two-dimensional 
Helmholtz  instability.  The  nonlinear  growth  of  this  instability 
leads  to  a rotational  two-dimensional  vortex  flow  with  the  vortices 
pairing  and  amalgamating  (see  Winant  and  Browand  and  Roshko). 
The  downstream  portion  of  Plate  la  and  b shows  the  flow  in  an 
amalgamation  state.  This  amalgamation  process  apparently 
produces  some  three-dimensionality.  This  three-dimensionality 
produced  by  pairing  can  be  seen  in  the  first  few  structures  shown 
in  Plate  lc,  d,  e and  f (and  Brown  and  Roshko  Fig.  20d  and  e). 
However,  high  speed  movies  of  this  process  and  pictures  of 


correspondingly  increasing  in  scale. 

However,  another  type  of  three-dimensionality  begins  to 
appear  just  upstream  of  the  critical  Reynolds  number  position. 
This  three-dimensionality  corresponds  in  the  photographs  to  the 


appearance  of  longitudinal  (streamwise)  lines  and  a "dimpling"  of 

the  surface  of  the  large  structures  (apparent  on  the  upstream  and 

downstream  edges  of  the  structure  in  the  plan-view  photographs). 

This  phenomenon  first  appears  at  a Reynolds  number  (&AU/V-.  ) 

4 

of  approximately  10  . A close  examination  of  the  plan -views 
suggests  that  the  three-dimensionality  is  cellular  in  nature  and 
that  these  structures  begin  as  hairpin  vortices  (Plate  If,  x = 1/3). 
Downstream  of  the  appearance  of  these  cellular  structures  as 
they  are  being  rolled  up  into  the  large  structures,  they  seem  to 
become  unstable  and  produce  the  three-dimensional  small  scales 
within  the  large  structures.  This  catastrophic  instability  of  the 
cellular  structures  always  first  takes  place  within  the  large 
structure  itself.  This  phenomenon  occurs  at  what  is  being  called 
the  critical  Reynolds  number.  Several  other  features  are  note- 
worthy. Once  these  cellular  structures  appear,  they  continue  to 
appear  in  all  downstream  structures  and,  in  fact,  seem  to  be 
continuous  through  several  structures  if  the  structures  are 
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approximately  equal  in  size.  Thus,  once  the  cellular  structures 
become  unstable  as  they  are  being  rolled  up  into  the  large  struc- 
tures, the  large  structures  become  characterized  by  a considerable 

18 

amount  of  fine  scale  mixing,  and  the  thin  vorticity  layers  (braids  ) 
connecting  the  large  structures  become  characterized  by  stream- 
wise  cellular  structures. 

At  even  larger  Reynolds  numbers  (Plate  lg  and  h),  the 
cellular  structure  instability  propagates  into  the  braids,  thus 
causing  much  of  the  layer  to  be  characterized  by  small  scale, 
three-dimensional  motions.  Even  though  there  is  an  increase  in 
the  amount  of  three-dimensionality  above  Re^,  the  shear  layer 
continues  to  grow  at  the  same  linear  rate  and  the  basically  two- 
dimensional  large  structures  continue  to  pair  and  grow.  The 
significance  of  this  will  be  explored  in  Section  6.  3. 

Before  proceeding  to  a possible  explanation  of  this  phenom- 
enon, the  results  of  a similar  photographic  investigation  of  other 
shear  flows  will  be  presented.  This  investigation  was  undertaken 
to  see  if  evidence  for  a similar  instability  in  the  shear  flows 
could  be  found. 

First  of  all,  in  order  to  verify  that  the  instability  is  not 
caused  by  the  density  difference  in  the  nitrogen-helium  shear 
layer,  the  constant  density  shear  layer  was  investigated.  Shown 
in  Plate  2 is  a Reynolds  number  survey  of  this  flow.  The 
shadowgraph  technique  is  still  usable  due  to  the  difference  in 
indices  of  refraction  between  nitrogen  and  the  mixture  of  argon 


and  helium. 
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The  similarities  between  this  case  and  the  last  are  striking. 
Again,  at  the  lower  Reynolds  numbers,  the  large  structures  are 
completely  two-dimensional  (Plate  2a  and  b,  up  to  x = 1/5  and 
Plate  2c  and  d,  up  to  x = 1/10).  At  a Reynolds  number  of  approx- 

4 

imately  0.  5 x 10  , the  cellular  structures  again  appear,  beginning 

4 

as  hairpin  vortices.  At  a Reynolds  number  of  about  1.5x10  , 
the  cellular  structures  become  unstable  within  the  large  structures 
producing  an  increase  in  fine  scale  mixing.  The  repetition  of  this 
process  is  clearly  shown  in  Plate  2e  and  f.  Here  the  cellular 
structures  are  quite  clearly  seen  in  the  braid  connecting  the  last 
two  structures,  both  of  which  show  fine  scale  mixing  produced  by 
the  instability  of  the  cellular  structures  within  the  large  structures. 
The  continued  basic  two-dimensional  nature  of  the  large  structures 
is  illustrated  in  Plate  2g  and  h.  The  spanwise  correlation  can 
clearly  be  seen  in  the  plan  view.  The  instability  of  the  cellular 
structures  within  the  braids  can  also  be  seen  throughout  much  of 
Plate  2h.  (The  apparent  free  stream  turbulence  is  due  to  slight 
cooling  of  the  gas  at  this  high  pressure.  ) Again,  it  can  be  seen 
that  the  growth  and  development  of  the  large  structures,  and  thus 
the  shear  layer  itself,  are  basically  unaltered  by  this  cellular 
structure  instability. 

The  next  case  studied  was  that  of  the  inhomogeneous  wake. 
This  case  is  quite  different  from  the  previous  cases  for  several 
reasons:  the  Reynolds  number  of  a wake  is  constant,  the  flow 

has  large  structures  of  opposite  vorticity,  the  similarity  growth 
laws  are  different,  and  it  is  not  clear  whether  or  not  a wake  has 
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an  amalgamation  process.  However,  in  spite  of  these  differences, 
the  inhomogeneous  wake  exhibited  a similar  instability.  One 
difference,  of  course,  between  the  instabilities  observed  is  that 
for  the  wake,  since  the  Reynolds  number  is  constant,  if  the 
cellular  structure  is  not  present  when  the  wake  develops,  it  will 
not  appear  later.  Likewise,  if  the  cellular  structure  is  present 
anywhere,  it  will  be  present  everywhere.  Plate  3a  and  b shows 
the  wake  at  a relatively  low  Reynolds  number.  The  completely 
two-dimensioiVal  nature  of  the  large  structures  is  evident.  Plate 
3c  and  d shoves  the  wake  at  a larger  Reynolds  number.  Again,  the 
basically  two-dimensional  nature  of  the  large  structures  is  seen; 
however,  the  three-dimensional  cellular  structure  is  now  present 
throughout  the  flow.  The  spanwise  correlation  is  still  evident  in 
the  slightly  higher  Reynolds  number  wake  shown  in  Plate  3e  and  f, 
but  now,  the  cellular  structure  has  become  somewhat  unstable. 

In  Plate  3g  and  h,  the  cellular  structure  has  become  quite  unstable. 
The  spanwise  correlation,  however,  is  still  remarkable.  Thus, 
an  instability  similar  to  the  one  observed  in  the  shear  layer  is 
present  in  the  inhomogeneous  wake. 

The  final  case  studied  was  a homogeneous  wake.  Plate  4a 
and  b shows  the  wake  before  the  cellular  structure  develops  and 
Plate  4c  and  d shows  the  wake  after  the  cellular  structure  has 
become  unstable.  Again,  a large  spanwise  correlation  is  seen  in 
both  cases  even  though  one  has  much  more  three-dimensionality. 

For  comparison,  the  momentum  thickness  was  greatly  increased 
for  the  final  photographs  of  this  case.  The  cellular  structure  is 


clearly  shown  in  the  upstream  portion  of  Plate  4e. 

6.  2 Explanation  of  Reynolds  Number  Effect 

To  review  what  has  been  observed  to  this  point,  the  concen- 
tration measurements  showed  that  above  a critical  Reynolds  number 
the  amount  of  small  scale  mixing  in  the  shear  layer  is  substantially 
increased.  The  entrainment  rates  into  the  turbulence  are  also 
larger  above  Re^,.  The  intermittency  profiles  showed  that  above 
this  Reynolds  number  the  probability  of  finding  pure  fluid  within 
the  layer  is  reduced.  Just  before  this  increase  in  small  scales 
occurs,  cellular  three-dimensional  structures  appear  in  certain 
portions  of  the  large  structures  beginning  as  hairpin  vortices. 

These  structures  become  unstable  as  they  are  being  rolled  up 
into  the  large  structures.  Thus,  the  large  structure  itself  has 
its  own  instability  and  above  a critical  Reynolds  number  this 
instability  generates  strong  local  three-dimensionality  and  a 
cascade  to  high  wave-numbers  within  the  basically  two-dimensional 
large  structure.  This  cellular  structure  instability  propagates 
into  the  braids  connecting  the  large  structures  thus  producing 
small  scale  three-dimensionality  throughout  much  of  the  layer. 

This  cellular  structure,  with  its  associated  instability,  repeats 
itself  at  all  Reynolds  numbers  above  Rer. 

Any  discussion  of  the  physical  characteristics  of  this 
instability  associated  with  the  large  structures  and  the  reasons 
for  its  continued  appearance  at  a Reynolds  number  above  the 
critical  value  is  somewhat  speculative  at  this  stage,  though  much 
can  be  inferred  from  the  photographs  and  concentration  data.  It 
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is  useful,  however,  to  construct  a model  which  is  at  least 

consistent  with  these  inferences,  if  only  to  provide  a possible 

basis  for  further  experiments.  In  addition,  since  the  development 

of  the  large  structure  is  essentially  an  unsteady  process  with  far 

from  simple  geometrical  features,  some  physical  appreciation  of 

the  mechanisms  involved  will  almost  certainly  be  required  before 

any  analytical  description  can  be  fruitfully  tackled. 

Since  this  new  structure  is  cellular  in  nature  and  appears 

to  be  associated  with  a secondary  flow,  it  is  more  than  likely 

that  this  instability  arises  basically  from  a Rayleigh  type  instability. 
1 9 

Rayleigh  showed  that  in  a steady,  inviscid,  purely  circulatory 
flow,  with  tangential  velocity  v(r),  the  radial  equilibrium  is 
unstable  if  the  square  of  the  circulation  decreases  outward. 

Since  ^ (r  v ) = 2r  uu  £ , where  C = 7 ^7  (rv)  is  vorticity  and 
uT  = v/r  is  angular  velocity,  this  inviscid  criterion  for  circulatory 
flow  can  also  be  expressed  in  the  form 

1 15  < 0 . 6.  1 

In  other  words,  the  motion  is  unstable  whenever  the  sense  of  the 

local  rotation  Q is  opposite  to  the  sense  of  the  general  rotation  u. 

Taylor^carried  the  analysis  a step  further  and  determined  the 

effect  of  viscosity  on  the  stability  boundary  both  experimentally  and 

theoretically. 

2 1 

Coles^  , using  dimensional  arguments  and  balancing  of  terms 


showed  that  the  viscous  Taylor  crite- 
2 

v k 7 
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where  fl  is  some  appropriate  characteristic  or  mean  value  for  the 
angular  velocity  uu(r)  and  L is  a characteristic  scale  for  the 
secondary  flow.  By  using  Taylor's  data,  he  showed  that  if  0 
is  approximated  by 

n = (ujj  + ^)/2  6,  3 

where  and  are  the  angular  velocities  on  the  boundaries  of 
the  region  of  instability,  then  L can  be  approximated  by 

i 

L = d/(1708)4  6.4 

where  d is  the  thickness  of  the  region. 

In  order  to  apply  this  criterion  to  the  present  case,  some 
estimation  of  the  terms  involved  must  be  made.  The  flow  field 
calculations  made  by  Corcos  and  Sherman  of  the  inviscid  roll-up 
of  a two-dimensional  shear  layer  can  be  used  to  make  these 
estimations.  Their  Figure  1 (reproduced  below)  shows  the 
streamline  pattern  and  region  of  vorticity  concentration  (shaded 
area)  for  this  roll-up  process  for  a reference  frame  moving  with 
the  structures. 

Since  the  regions  marked  B and  C have  vorticity  of  opposite 
sense  to  the  general  rotation  within  the  regions,  it  is  apparent 
that  these  regions  meet  the  Rayleigh  criterion  (Eq.  6.1)  for 
instability.  Therefore,  rough  estimations  of  the  terms  in  Eq.  6.  2 
will  be  made  for  these  regions  to  determine  the  conditions  required 
for  these  regions  to  also  be  unstable  according  to  the  Taylor 
criterion.  All  of  the  estimations  to  be  made  are  quite  rough  and 
only  made  to  see  if  a Taylor-type  instability  is  consistent  with 
the  observed  conditions  within  the  large  structures.  What  is 


SENSE  OF  VORTICITY-O 


SECTION  A-A 


suggested  here  is  that  the  Taylor  instability  gives  rise  to  hairpin 
vortices  in  which  the  "head"  of  the  hairpin  originates  in  the 
unstable  region  of  a large  structure  and  is  connected  to  longitu- 
dinal vortex  lines  in  the  braid.  Such  a flow  structure  is  shown 
schematically  as  Section  A-A  in  the  drawing.  Two  pairs  of  Taylor 
vortices  are  shown.  Such  a flow  field  would  show  up  in  the  plan- 
view  as  several  parallel  lines  (four  are  shown)  running  in  the 


46 


streamwise  direction  in  the  braid.  These  lines  would  be  produced 

by  the  interfaces  between  the  two  gases.  If  d is  the  thickness  of  the 

braid  connecting  the  large  structures,  then  the  mean  vorticity  (Q 

can  be  approximated  by  AU/d.  Since  the  center  of  the  braid  is  a 

stagnation  point,  can  be  taken  to  be  zero.  Therefore,  if  the 

radius  of  curvature  of  the  streamline  on  the  edge  of  region  B is 

taken  to  be  lOd,  then  can  be  approximated  as  -(AU/2)/10d. 

Thus,  according  to  Eq.  6.3,  Q l -(A  U/2)/l  0 d/2  = -AU/40d.  Since 

the  width  of  the  region  is  d/2,  L becomes  (d/2)/(1708)4.  Upon 

inserting  these  approximations,  Eq.  6.2  becomes 

2 


M)(M)  > 
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\( 1708)*/ 


^ > 740  . 6.5 

v 

In  Taylor's  experiment  it  was  found  that  the  spacing 
between  vortices  (dj)  is  approximately  equal  to  the  region  thickness 
(d).  From  the  photographs  taken  in  this  investigation,  d^  was 
found  to  be  approximately  one -tenth  of  the  visual  thickness  of  the 
layer  (^y^g)*  Since  6^  is  approximately  one -half  of  6y^s*  d/6^ 
is  approximately  0.  2.  Therefore,  the  Taylor  criterion  for 


instability  becomes 


6 AU  . 

— — > 0.  4 x 104 

v 


The  closeness  of  this  estimation  of  the  point  of  onset  of  the 
Taylor  instability  to  the  observed  point  of  onset  is  more  coinci- 
dental than  remarkable.  Other  equally  valid  estimations  of  the 
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terms  involved  would  have  obviously  given  different  answers. 
However,  what  this  dees  demonstrate  is  that  the  characteristics 
o1  the  instability  observed  in  this  experiment  are  remarkably 
similar  to  the  characteristics  of  the  Taylor  instability  and  that  the 
conditions  required  for  a Taylor  instability  are  present  in  some 
regions  of  the  large  structures.  It  is  also  significant  that  the 
onset  of  the  Taylor  instability  occurs  at  a lower  Reynolds  number 
in  the  homogeneous  case  than  in  the  inhomogeneous  case.  The 
kinematic  viscosity  used  for  nondimensionalization  in  the  inhomo- 
geneous shear  layer  should  actually  be  an  average  in  some  way  of 
arid  v^e  rather  than  just  . Thus,  since  any  average  of 
Vj^  and  v^e  will  be  larger  than  , a Reynolds  number  based  on 

V..  will  be  larger  than  a Reynolds  number  for  the  same  conditions 
based  on  an  average  v.  Thus,  the  Reynolds  number  of  the  onset 
of  the  Taylor  instability  will  appear  to  be  higher  if  the  Reynolds 
number  is  based  on  rather  than  an  average  v.  Plate  If  further 

demonstrates  this  kinematic  viscosity  dependence.  The  Taylor 
instability  first  appears  on  the  nitrogen  side  of  the  large  structure 
(region  B)  rather  than  on  the  helium  side.  This  behavior  is  in 
accord  with  Eq.  6.  6. 

The  large  structure  not  only  provides  the  conditions 
necessary  for  the  development  of  this  instability,  but  it  also 
provides  a mechanism  to  amplify  it.  Once  the  longitudinal  vorticity 
produced  by  the  Taylor  instability  is  present  in  the  braid,  the 
vortex  stretching  due  to  the  high  strain  rate  in  the  plane  of  the 
braid  (as  shown  in  the  drawing)  strengthens  this  vorticity. 
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Apparently,  as  these  Taylor  vortex  pairs,  which  have  been 
strengthened  by  the  stretching  in  the  braid,  are  being  rolled  up 
into  the  large  structure,  the  additional  vortex  stretching  causes 
the  pairs  to  precess  around  each  other  and  eventually  break  up 
into  small  scale  structures.  As  can  be  seen  in  the  photographs, 
the  generation  of  small  scales  is  most  apparent  on  the  leading 
convex  surface  of  the  structure,  on  the  high  velocity  side  of  the 
layer,  and  on  the  corresponding  trailing  surface  on  the  low 
velocity  side,  which,  from  the  above  picture,  is  where  the  Taylor 
vortices  are  expected  to  break  down.  The  precise  mechanism  by 
which  these  secondary  flows  break  down  to  small  scale  motions 
is  not  entirely  clear,  but  what  is  clear  is  that  this  motion  is 
unstable  in  the  sense  that  the  vortex  stretching  in  the  braids  and 
the  Taylor  instability  in  the  large  structures  are  coupled  so  that 
there  is  a strong  mechanism  which  can  transfer  energy  from  the 
mean  flow  into  three-dimensional  motions. 

Two  final  points  concerning  the  spacing  between  pairs  of 
Taylor  vortices  (d^)  should  be  noted.  First  of  all,  the  downstream 
disturbance  produced  by  the  Taylor  vortices  in  one  braid  is 
apparently  sufficient  to  determine  the  positions  of  the  vortices  in 
the  next  braid.  Since  both  braids  are  connected  (in  a sense) 
through  one  large  structure,  it  is  quite  easy  to  see  how  this 
communication  of  position  can  take  place.  This  upstream 
determination  of  the  position  of  the  vortices  is  seen  in  the  photo- 
graphs as  the  apparent  continuation  of  the  vortices  through  several 
large  structures.  The  second  point  concerns  the  absolute  value 
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of  the  spacing  between  the  pairs  of  vortices.  In  the  mean, 
was  found  to  be  approximately  one-half  6 . However,  the 

downstream  influence  just  discussed  is  sufficiently  strong  so  that 
d^  changes  in  a stepwise  manner.  Thus,  in  some  cases,  d^ 
is  such  that  40,  20,  or  10  pairs  can  fit  across  the  width  of  the 
layer  (10  cm),  and  in  other  cases,  28,  14,  or  7 pairs  can  fit. 

In  other  words,  the  spacing  is  constant  until  the  point  is  reached 
at  which  approximately  half  the  number  of  pairs  across  the  layer 
is  more  stable.  It  is  therefore  apparent  that  the  spacings  between 
pairs  of  vortices  is  a function  both  of  6 and  the  spanwise  boundary 
conditions. 

6.  3 Implications  of  the  Reynolds  Number  Effect 

Asa  final  comment  on  this  three-dimensional  instability, 
some  of  the  important  implications  of  this  phenomenon  will  be 
pointed  out. 

As  was  discussed  in  the  introduction  to  this  work,  the 
basically  two-dimensional  large  structure  observed  in  the  shear 
layer  was  thought  to  not  be  able  to  undergo  vortex  stretching  and 
therefore  could  not  form  a part  of  the  fully-developed  energy 
cascade  from  mean  flow  to  smallest  eddies,  an  essential  character- 
istic of  turbulence.  Now,  with  the  observed  combination  of 
production  of  Taylor  vortices  and  vortex  stretching,  the  process 
by  which  the  large  structure  becomes  a part  (in  fact,  an 
indispensable  part)  of  the  cascade  to  small  scales  is  clearly  seen. 
This  process  is  continuously  repeated  in  the  development  of  all 
large  structures  above  the  critical  Reynolds  number  and  is  a 
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basic  mechanism  in  the  continued  generation  of  three-dimensional 
turbulence  energy.  As  evidenced  by  the  pictures,  there  is  little 
indication  of  a cascade  to  small  scales  below  the  critical  Reynolds 
number.  The  high-wave-number  energy  in  the  power  spectral  den- 
sity curves  below  Re^.  comes  from  the  sharp  interfaces  in  the  flow, 
not  from  small-scale  spatial  structures. 

A further  consequence  of  this  phenomenon  can  be  seen  by 
what  the  three-dimensional  instability  does  not  do;  it  does  not 
significantly  affect  the  linear  growth  rate  of  the  shear  layer.  This 
of  itself  is  quite  significant  and  demonstrates  the  fundamental  role 
which  the  large  structure  plays.  In  the  flow  both  above  and  below 
the  critical  Reynolds  number,  the  Reynolds  stress  clearly  derives 
from  the  large  structure.  However,  above  the  critical  Reynolds 
number,  a further  mechanism  for  the  production  of  turbulence 
exists  which  is  not  at  all  related  to  the  Reynolds  stress.  Thus, 
the  Reynolds  stress  is  not  a complete  measure  of  the  production 
of  turbulence  energy. 

As  remarked  earlier,  the  entrainment  rates  into  the  turbu- 
lence are  increased  by  this  instability  phenomenon.  Due  to  the 
breakdown  of  the  Taylor  vortices,  the  active,  newly  formed  turbu- 
lence on  the  leading  and  trailing  convex  surfaces  of  the  large 
structure  evidently  enhances  the  local  entrainment  rates  at  these 
surfaces.  This  reduces  the  penetration  of  irrotational  fluid  into 
the  layer  as  can  be  seen  from  the  intermittency  plots.  Thus, 
above  the  critical  Reynolds  number,  the  entrainment  of  the  fluid 
into  the  turbulence  is  partially  by  "engulfment"  due  to  the  induced 


velocity  field  of  the  large  structure  and  the  amalgamation  process 
and  partially  by  nibbling.  The  fact  that  the  entrainment  rate  into 
the  turbulence  increased  by  only  8%  above  the  critical  Reynolds 
number  suggests  that  nibbling  is  of  lesser  importance  in  the  flow 
than  engulfment. 

The  existence  of  these  secondary  instabilities  helps  to 
account  for  the  large  v'wr  correlations  which  have  been  observed 
by  some  investigators  . A completely  two-dimensional  view  of 
these  large  structures  cannot  account  for  these  large  correlations. 

Finally,  it  is  noted  that  all  of  these  conclusions  have  been 
drawn  from  experiments  on  a mixing  layer  having  a velocity  ratio 
of  r = U-/7 . There  is  no  reason,  however,  to  doubt  that  the 
results  are  applicable  to  mixing  layers  of  any  velocity  ratio.  The 
existence  of  large  coherent  structures  in  wakes  at  very  high 
Reynolds  numbers  and  the  strong  evidence  of  a similar  instability 
in  the  homogeneous  and  inhomogeneous  wakes  shown  in  the  photo- 
graphs of  these  flows  suggests  that  this  instability  may  be  an 

integral  part  of  the  production  of  turbulence  energy  in  many  other 

* 

shear  flows. 


*Many  of  the  ideas  presented  in  this  Section  concerning  the  explana- 
tion of  the  observed  Reynolds  number  effect  were  developed  in  dis- 
cussions with  Dr.  Garry  Brown  and  Professor  Anatol  Roshko.  1 
am  especially  indebted  to  Dr.  Brown  for  the  many  hours  he  spent 
thinking  about  this  phenomenon  and  sharing  his  ideas  with  me. 
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VII.  APPLICATIONS  TO  CHEMICALLY  REACTING  SHEAR  FLOWS 

i 

7.  1 Discussion  of  Toor's  Analysis 

The  unmixedness  results  discussed  in  Section  5 are,  for 

the  most  part,  qualitative.  Little  more  can  be  said  in  applying 

these  results  to  chemically  reacting  shear  flows  than  that  one 

would  expect  more  product  produced  for  one  flow  than  for  the 

other.  However,  under  certain  restrictions,  the  concentration  data 

obtained  can  be  used  to  predict  reaction  rates  for  some  special 

cases.  An  analysis  of  reactions  under  these  conditions  was  first 
23 

put  forth  by  Toor  . His  analy  is  shows  that  if  a few  assumptions 
are  made,  the  probability  density  functions  of  the  concentrations  of 
scalar  contaminants  in  a nonreacting  flow  can  be  used  to  predict 
the  amount  of  product  formed  in  a reacting  flow  of  similar  geometry. 

This  restriction  of  similarity  between  the  reacting  and  nonreacting 
flows  requires  that  the  reaction  be  dynamically  passive.  This  will 
be  true  if  the  reactants  are  dilute  and  the  reaction  releases  only 
small  amounts  of  heat. 

The  reacting  and  nonreacting  flows  for  this  analysis  are 

illustrated  in  Figure  1.  CLq  is  the  free- stream  concentration  of 

each  reactant  (CLq  « 1).  In  the  nonreacting  flow,  the  concentration 

of  fluid  1 (Cu)  is  measured  as  a function  of  position  and  time.  In 
ri 

the  reacting  flow,  reactant  A is  assumed  to  be  diluted  in  fluid  1 
and  reactant  B in  fluid  2;  A and  B are  assumed  to  react  irrevers- 
ibly by  the  reaction 

NaA  + NgB  - NpP  . 7.  1 

A 
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In  order  for  this  analysis  to  be  valid,  the  diffusivities  of  all 
species  in  both  flows  must  be  equal  (D^  = Dg  = Dp  = Dj  = D£). 
Also,  the  reaction  rate  must  be  much  faster  than  the  diffusion 
rate.  Many  of  these  restrictions  can  be  considerably  relaxed 
before  large  errors  are  introduced. 

The  basic  idea  behind  Toor's  analysis  as  applied  to  the  shear 
layer  can  be  stated  quite  simply.  Once  the  relative  concentration 
of  the  two  fluids  in  the  nonreacting  case  are  known  for  some  posi- 
tion, then,  given  the  free- stream  concentrations  of  reactants  for 
the  reacting  case,  the  relative  concentrations  of  the  reactants  can 
be  found  for  that  position.  However,  since  fast  chemistry  is  as- 
sumed, the  two  reactants  can  not  co-exist;  they  must  react.  Thus, 
for  a given  stoichiometry,  the  reactant  that  is  in  excess  can  be 
determined;  the  other  reactant  concentration  will  be  equal  to  zero 
and  the  concentrations  of  the  product  and  of  the  excess  reactant  can 
be  calculated.  Thus,  the  P.  D.  F.  's  of  the  reactant  and  product  con- 
centrations can  be  found  by  calculating  these  concentrations  for  each 
scalar  concentration  and  using  the  scalar  P.  D.  F.'s  to  determine  what 
percentage  of  time  each  occurs.  The  complete  methodology  of  applying 
Toor's  analysis  to  this  investigation  is  given  in  Appendix  D.  It  was 
found  to  be  convenient  to  normalize  all  concentrations  with  the  corres- 
ponding stoichiometric  coefficients.  Thus, 

= (molar  concentration  of  species  i)/N..  7.2 

7.2  "Equivalent  Laminar  Solution" 

Another  often  used  method  for  predicting  the  amount  of 
product  in  a reacting  shear  flow  is  what  can  be  called  an 


"equivalent  laminar  solution"  . This  method  uses  only  mean 
profiles  and  neglects  all  concentration  fluctuations.  It  is  as- 
sumed that  all  fluid  which  is  entrained  into  the  layer  is 
immediately  mixed  on  a molecular  level.  Thus,  for  a given 
entrainment  rate  and  a mean  concentration  profile,  it  is  an 
upper  bound  on  the  amount  of  product  which  can  be  produced. 
By  this  analysis,  (the  concentration  of  reactant  A as  cal- 

culated by  the  laminar  solution)  would  be  equal  to  CAoCH  and 

C„T  would  be  equal  to  Cn  (1  - C,,)  except  for  the  fact  that 
B L 130  n 

the  two  reactants  can  not  co-exist.  Thus,  one  concentration 
must  be  reduced  by  the  reaction  and  the  other  completely  re- 
acted away.  Therefore, 


c . „ 

= C . C..  - 
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= CBo  <>  - CH>  lf  CAL  > 0 ' 

These  equivalent  laminar  solutions  will  be  used  mainly  to 
show  the  effect  of  the  unmixedness  of  the  gases  on  the  amount 
of  product  formed. 
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7.  3 Results  and  Discussion 

Discussed  in  this  section  are  the  results  of  applying  the 
analysis  of  Sections  7.1  and  7.2  to  the  three  cases  investigated. 

Figure  18  shows  the  results  of  these  calculations  for  the 
inhomogeneous  shear  layer  at  a Reynolds  number  above  Re^.  The 
mean  concentration  profiles  are  shown  in  Figures  18a  and  c and 
the  fluctuation  profiles  are  shown  in  Figures  18b  and  d.  Two  dif- 
ferent cases  are  shown  here  for  comparison.  In  one  case  the 
ratio  (a)  of  the  entrained  reactant  A to  the  entrained  reactant  B is 
set  equal  to  1,  and  in  the  other  case,  this  ratio  is  set  equal  to  3. 
The  actual  definition  of  a is 


a = 


Ve, 


Ve. 


7.6 


Thus,  the  effect  of  changing  freestream  concentrations,  entrain- 
ment rates,  or  stoichiometric  coefficients  (see  Eq.  7.  2)  can  be 
shown  simply  by  changing  a. 

The  effect  of  varying  a can  be  seen  by  comparing  Figure  18a 
with  Figure  18c.  For  a = 1,  since  Ve^  = 3\Te^,  CgQ  must  be 
about  three  times  as  great  as  Thus,  since  CgQ  was  unity 

for  all  calculations,  CA/Cgo  approaches  approximately  0.3  on  the 
helium  side  and  Cg/CgQ  approaches  1.  0.  Obviously,  since  there 
is  more  reactant  available  for  a = 3 than  for  a = 1,  Figure  18c 
shows  a higher  concentration  of  product. 

The  equivalent  laminar  solution  profiles  are  shown  for  the 
two  cases  by  the  dashed  curves.  The  difference  between  the 


laminar  solution  and  the  P.  D.  F.  solution  shows  the  effect  of  the 
unmixedness  of  the  gases.  For  the  laminar  solution,  A and  B 
cannot  exist  at  the  same  r\  location.  Therefore,  the  profiles  go  to 
zero  at  the  same  location. 

The  RMS  concentration  fluctuations  are  shown  in  Figures  18c 
and  d.  is  normalized  with  the  value  of  Cp  calculated  for 

the  laminar  solution  for  the  point  at  which  both  CA  and  CR  go  to 
zero. 

The  results  for  the  inhomogeneous  shear  layer  at  a Reynolds 
number  less  than  Re^  are  shown  in  Figure  19.  Again,  two  cases 
are  shown:  a = 1 and  a = 3.  The  major  differences  to  note  be- 

tween Figures  18  and  19  are  that  C^,  is  lower  and  the  fluctuations 
are  higher  in  Figure  19  than  in  Figure  18. 

The  homogeneous  shear  layer  results  are  shown  in  Figure 
20,  and  those  for  the  wake  are  shown  in  Figure  21. 

The  effect  that  Reynolds  number  has  on  the  extent  of 
molecular  mixing  is  shown  in  Figure  22.  Plotted  in  Figure  22a  is  the 
integral  of  the  product  concentration  normalized  with  the  integral 
of  the  laminar  product  concentration  as  a function  of  Reynolds 
number  for  the  inhomogeneous  shear  layer  P a*  1 P^  are  defined 


The  rather  sharp  increase  in  the  amount  of  product  at  the  critical 


mu 
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Reynolds  number  again  illustrates  the  increase  in  the  extent  of  molecu- 
lar mixing  above  Rp^,.  It  is  evident  from  this  plot  that  the  extent  of 
molecular  mixing  is  constant  both  above  and  below  the  critical  Reynolds 
number,  the  mixing  above  Rec  being  approximately  25%  greater  than 
that  below  lie-.  P/PT  for  the  homogeneous  shear  layer  was  found  to 
be  0.  75  and  for  the  wake,  0.83.  These  values  are  for  conditions  above 
the  critical  Reynolds  number  estimated  from  the  shadowgraphs  for 
each  case. 

Since  P has  the  dimensions  of  a length,  it  can  be  thought  of  as 
a measure  of  the  "width"  of  the  shear  layer.  Therefore,  to  show  the 
effect  that  Reynolds  number  has  on  the  growth  rate  based  on  this 
width,  the  Reynolds  number  based  on  P can  be  plotted  as  a function  of 
the  Reynolds  number  based  on  5.  This  is  shown  in  Figure  22b.  The 
laminar  solution  is  shown  for  comparison. 

Some  insight  into  a rather  interesting  problem  can  be  gained  by 
a further  application  of  the  results  of  this  analysis.  The  problem  is  to 
determine  how  much  excess  of  one  reactant  is  required  in  order  to 
completely  react  all  of  the  other.  For  example,  in  this  case,  suppose 
reactant  A is  an  oxidant  and  reactant  B is  a fuel.  Then,  the 
problem  is  to  determine  what  freestream  concentration  of  oxidant 
is  required  so  that  all  of  the  entrained  fuel  will  react.  Asa  larger  and 
larger  percentage  of  fuel  is  being  reacted,  the  variation  of  the  flux  of 
product  past  any  streamwise  station  will  approach  the  value  of  entrained 
fuel  at  that  station.  Therefore,  the  required  excess  of  oxidant  can  be 
determined  by  finding  the  ratio  of  the  streamwise  variation  of 
product  flux  to  the  entrained  fuel  as  a function  of  the  ratio  of 
entrained  oxidant  to  entrained  fuel.  Thus,  if  the  flux  of  product 
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is  defined  as 


d ~ 

what  is  required  is  a plot  of  CgQ^  (PJ/fV^C^/Uj ) as  a function 

of  (Ve,CA  /Ve_C„  ).  The  results  of  this  analysis  are  shown  inFig- 

ure  23.  Again,  the  Reynolds  number  effect  is  seen  quite  clearly 

in  Figure  23a  which  shows  the  results  for  the  inhomogeneous  shear 

layer  both  above  and  below  the  critical  Reynolds  number.  The 

results  for  the  homogeneous  shear  layer  and  the  wake  are  shown 

in  Figures  23b  and  2 3c,  respectively.  Figure  23d  shows  all  four 

curves  together  for  comparison. 

The  basic  result  of  this  analysis  is  that  approximately  ten 

times  as  much  oxidant  must  be  entrained  as  fuel  in  order  to  completely 

react  the  fuel.  Thus,  for  the  inhomogeneous  shear  layer  and  for  unit 

stoichiometric  coefficients,  C.  would  have  to  be  a little  over 

o 

three  times  as  great  as  CgQ(since  Ve^/Ve^  = 3).  Likewise,  for  the 

wake,  C.  must  be  approximately  seven  times  as  great  as 
o 

Cgo(since  Ve^/V^  = 3/2). 

From  the  curves  in  Figure  23,  the  streamwise  dependence 
of  P can  be  determined.  In  other  words, 

(fJ)SL  = bSL(x  - Xo>  7'9 

for  the  shear  layer  and 

(l).  ■ 

for  the  wake  can  be  found,  b.  is  of  course  a function  of  a. 


Thus,  knowing  this  information,  a comparison  between  the  wake 


and  the  shear  layer  can  be  made.  Such  a comparison  will 


determine  which  case  produces  more  flux  of  product  under  similar 


conditions.  Obviously,  the  answer  depends  on  the  x station  at  which 


the  comparison  is  being  made.  As  discussed  in  Section  5.  3c,  it 


is  expected  that  the  wake  will  initially  have  a larger  product  flux 


since  it  has  better  molecular  mixing  and,  initially,  more  entrain- 


ment than  the  shear  layer.  This  analysis  will  determine  the  x 


station  at  which  the  shear  layer  begins  to  have  a larger  product 


flux.  In  order  for  the  comparison  to  be  equitable,  the  reactant 


fluxes  were  arranged  to  be  the  same  for  both  cases.  This 


required  that  (U]CAo)SL  = <UiCAo*W  and  *U2CBo*SL  " ^U2CBo)W 
Thus,  for  (Uj)gk=  the  freestream  concentrations  would  be 


<CAo>SL  = <CAo»W  = '•»  »»<•  «CBo»SL  = <CBo>W/0‘ 38  = 


Therefore,  since  (Ve^/V^^^  = ^ anc*  ^Vei^e2^W  = °SL  = 


3. 4 and 


•=  3.  7.  For  these  values  of  a 


CBoii|P»SL/|Ve2CBo/Ul»  -°- 


CBo^(P>W/<Ve2CBo/Ul>  -0'11 


Therefore, 


CBo(P)SL  = °‘018<x-V 


• i. 

g /X  - X \ 2 

cbo<I»w  =°-091(-r-°)  • 


Finally,  the  x station  at  which  the  product  flux  will  be  equal  for 


both  cases  will  be 


256 


7.  11 
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or,  for  this  case,  (x  - x ) = 2.  8 inches.  This  means  that  for 
these  freestream  concentrations  of  reactants,  CgQ(P).yy  will  be 
larger  than  C (P)eT  up  to  (x  - x ) = 2.  8 inches. 

jdO  oL  O 

For  the  range  of  a investigated  (0.  1 < a<  10),  the  following 
empirical  relationships  between  P and  P were  found: 
for  the  inhomogeneous  shear  layer  (Re^.+  ), 

(P/P)  = 1.  73  + 0.  13  log  a , 7.  12 

for  the  inhomogeneous  shear  layer  (Re^."), 

(?/?*)  = 1.  76  + 0.  13  log  a , 7.13 

for  the  homogeneous  shear  layer, 

(P/P)  = 1.  39  + 0.  11  log  a,  7.  14 

and  for  the  wake, 

(P/P)  = 1.  085  - 0.  015  log  a.  7.  15 

Therefore,  by  combining  the  information  in  Figure  23  with  that  in 
Equations  7.  12  and  7.  15,  this  same  comparison  between  the  shear 
layer  and  the  wake  can  be  made  on  the  basis  of  P instead  of  P. 
When  this  was  done  it  was  found  that  CgQP  is  larger  for  the  wake 
than  for  the  shear  layer  up  to  (x  - xq)  = 99  (1  inch). 

One  final  comparison  was  made  in  which  all  freestream 
concentrations  were  set  equal  ((cAo)SLi=  ^Ao^^Bo'sL'^Bo'w1, 

Ami 

It  was  found  from  this  analysis  that  C^P  Is  larger  for  the  wake 
than  for  the  shear  layer  up  to  (x  - xq)  = 10.9  inches  and  CBqP  is 
larger  up  to  (x  - xq)  =3.9  inches.  This  increase  in  the  trade 
off  x station  is  of  course  due  to  the  increase  in  freestream  concen- 
tration flux. 

Curves  similar  to  those  in  Figure  2’  were  suggested  by  Garry 


Brown. 
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VIII.  CONCLUDING  REMARKS 

This  investigation  was  able  to  meet,  to  some  extent,  the 
goals  that  were  stated  in  the  Introduction  to  this  work.  The 
extent  of  molecular  mixing  was  measured  for  several  turbulent 
shear  flows.  It  was  determined  that  much  of  the  engulfed  fluid 
is  still  in  a high  state  of  unmixedness.  Even  the  fluid  which  has 
been  entrained  into  the  turbulent  regions  is  far  from  being  perfectly 
mixed  on  the  molecular  level.  It  was  further  found  that  above  a 
certain  critical  Reynolds  number,  the  extent  of  molecular  mixing 
increased.  Accordingly,  the  fluctuations  in  the  concentration 
decreased.  These  results  were  made  more  quantitative  by  using 
the  data  to  predict  the  results  when  a diffusion-limited  chemical 
reaction  is  incorporated  into  the  flow.  The  results  of  this  analysis 
showed  the  effect  that  the  lack  of  molecular  mixing  has  on  the 
progress  of  a chemical  reaction. 

Further  understanding  of  the  growth  and  development  of 
the  large  structures  in  turbulent  shear  flows  came  from  a rather 
unexpected  source.  A three-dimensional  instability  associated 


with  the  large  structure  itself  was  discovered  while  attempting 


to  find  the  cause  of  the  effect  varying  the  Reynolds  number  has 


on  the  extent  of  molecular  mixing.  It  may  now  be  seen  that  the 
large  coherent  structures  are  not  only  an  integral  part  of  the 
mechanism  which  controls  the  growth  of  the  shear  layer  but  also 
play  a central  role  in  the  production  of  small  scale,  three-dimen- 
sional turbulence. 
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APPENDIX  A 


Concentration  Probe 


The  probe  being  used  to  measure  concentration  histories 

in  turbulent  mixing  regions  is  a development  of  the  Brown-Rebollo 

probe  in  which  the  resolution  has  been  improved  by  reduction  in 

_5 

size.  This  is  accomplished  by  stretching  a hot-wire  (5  x 10 
inches  in  diameter)  across  the  diameter  of  a small  (0.  2 mm) 
capillary  tube  (see  Fig.  2a).  The  capillary  tube  is  necked  down 
upstream  of  the  wire  to  a very  small  throat  (0. 01  mm)  so  that 
when  a vacuum  is  applied  to  the  downstream  end  of  the  capillary 
the  flow  in  the  throat  is  sonically  choked.  The  flow  past  the  wire 
is  then  independent  of  the  freestream  velocity  (to  order  mJ*)  and 
is  dependent  only  on  the  throat  conditions  (y*,  a*,  pU*).  The 
conditions  at  the  wire  are  thus  functions  only  of  the  properties  of 
the  gas  being  sampled. 

The  probe  is  calibrated  by  placing  it  in  a chamber  in 

which  the  pressure  (P  ) and  the  relative  concentrations  (C)  of  the 

o 

two  gases  to  be  used  in  the  experiment  can  be  controlled.  Thus, 
the  voltage  output  ^proije  = V(PQ,  C)  is  determined.  These  data 
are  then  curve  fitted  to  give  C = C(Vp,  PQ). 

Estimates  of  the  spatial  resolution  of  the  probe  can  be 
made  by  calculating  the  captured  streamtube  diameter,  which  is 
a function  of  the  sonic  throat  diameter  and  the  freestream  Mach 
number.  For  example,  for  a freestream  velocity  of  1000  cm/'sec 
and  a throat  diameter  of  0.  01  mm,  the  streamtube  diameter  (and 
the  lateral  resolution)  is  approximately  0.  1 mm.  The  streamwise 


spatial  resolution  is  a function  of  the  hot-wire  frequency  response 
and  the  freestream  velocity.  Again,  for  a freestream  velocity 
of  1000  cm/sec  and  a wire  response  of  100  kHz,  the  streamwise 
spatial  resolution  is  0.  1 mm.  The  wire  response  is  estimated 
from  a standard  square  wave  test  for  a constant  temperature 
hot-wire.  A sample  square  wave  test  response  is  shown  in  Fig. 
2b.  An  interpretation  of  this  response  signal  can  be  found  in 
Ref.  25. 

An  operational  calibration  of  the  probe  response  was 
obtained  by  making  use  of  the  concentration  signal  (C(t))  itself. 
Interfaces  in  the  mixing  region  provide  very  high  values  of  dC/dt 
at  the  probe.  Analysis  of  the  frequency  spectra  of  the  signals 
obtained  under  various  conditions  (6,  U and  x can  be  varied)  shows 
whether  any  differences  can  be  attributed  to  resolution  limitations 
of  the  probe.  In  this  way  it  was  concluded  that  the  probe  has  a 
frequency  response  of  at  least  45  kHz,  which  corresponds  to  a 
spatial  resolution  of  0.2  mm  at  Uj  = 1000  cm/sec. 


*Much  of  the  work  done  to  improve  the  concentration  probe  was  done 
by  Dr.  Garry  Brown. 
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APPENDIX  B 

Minimization  of  the  Effect  of  Hydrostatic  Forces 
In  the  inhomogeneous  shear  flows  studied,  the  difference 
in  the  densities  of  the  two  fluids  produces  hydrostatic  forces 
which  induce  streamwise  gradients  in  both  pressure  and  velocity. 
The  magnitude  of  this  effect  can  be  seen  by  a simple  application 
of  the  steady  Bernoulli  equation,  where  ( is  air, 

= Pc  + |PU2  + (PA-P)gx  . B.  1 

OS 

The  last  term  is  included  since  the  flow  is  vertical  and,  therefore, 
the  relative  hydrostatic  forces  must  be  taken  into  account  when 
measuring  pressures.  Thus,  the  coefficient  of  pressure  for  the 
difference  in  pressure  between  two  vertical  points  (subscripts  a 


and  b)  is 


ps  - Ps 
b a 


$- 


If  the  fluid  is  helium,  then 


1 - H 


and  if  the  fluid  is  nitrogen,  then 

C = 1 - N2  - S , B, 

N 

where  H is  the  ratio  of  the  velocities  at  the  two  points  (U^/I^) 
for  helium,  N is  (U^/U  ) for  nitrogen  and 


Since  this  experiment  consists  basically  of  two  adjacent  columns 


! | 

I ; 


’Vs’Vf? 


of  gas  flowing  in  contact  with  each  other,  the  static  pressure  will 


adjust  itself  (through  changing  velocities)  so  that  the  pressure  will 

be  approximately  uniform  at  any  streamwise  station.  Therefore, 

CDtt  will  be  equal  to  since  p^,U,  = pT  U for  the  inhomo - 

rN  N H H 

geneous  shear  layer.  Thus,  Eqs.  B.  3 and  B.  4 can  be  combined 


to  give 


1 - H2  - SH  = 1 - N2  - SN 


N2  - H2  = SH  - SN  = S 


Equations  B.  3,  B.  4 and  B.  6 show  to  what  degree  the  hydrostatic 

forces  will  cause  gradients  in  the  pressure  and  velocity.  For 

example,  if  C and  C_  were  adjusted  to  be  zero  (by  positioning 
PH  PN 

the  side  walls)  and  if  the  pressures  measured  were  for  ports  at 
streamwise  stations  10  cm  apart,  then  the  velocities  could  change 
by  a total  of  more  than  5%  for  low  speed  flows. 

Since  constant  freestream  velocities  are  desired,  the  static 
pressure  should  be  adjusted  so  that  N and  H are  as  close  to  unity 
as  possible.  In  other  words,  the  pressure  gradient  should  be 
adjusted  to  minimize  the  sum 

(1  - H)2  + (1  - N)2  . B.  7 

Upon  eliminating  H in  Eq.  B.  7 by  using  Eq.  B.  6 and  using  the 
standard  minimization  technique,  the  following  fourth  order 
equation  for  N is  found: 

N4  - N3  + N2S  - NS+  S/4  = 0 . B.  8 


The  proper  root  of  Eq.  B.  8 can  then  be  substituted  into  Eq.  B.  4 
to  determine  what  pressure  gradient  will  minimize  the  effect  that 
the  hydrostatic  forces  have  on  the  velocities.  By  using  this 


t 


technique,  the  decrease  in  freestream  velocity  on  the  helium  side 
could  be  kept  to  less  than  2#  for  Ua  = 1000  cm/ sec  and  less  than 
0.5^  for  U = 2000  cm/sec.  The  corresponding  nitrogen  free- 
stream  velocities  increased  by  less  than  2$  and  0.  5$.  A similar 
procedure  was  used  for  the  inhomogeneous  wake. 


ifc-''  '-I  .-irli-  i 
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APPENDIX  C 

Technique  Used  to  Measure  the  Density  of  a Mixture  of  Gases 

In  order  to  study  the  mixing  process  in  a constant  density 
shear  layer,  two  gases  of  equal  density  but  different  concentration 
probe  response  were  needed.  Two  such  gases  are  nitrogen  and 
a homogeneous  mixture  of  1/3  helium  and  2/3  argon.  Since  the 
mixing  of  the  helium  and  argon  could  not  easily  be  controlled, 
some  means  of  measuring  the  density  of  the  mixture  was  needed 
in  order  to  determine  if  the  mixture  did  indeed  have  the  same 
density  as  nitrogen. 

Such  a measurement  can  simply  be  made  by  comparing  the 
hydrostatic  pressure  of  the  mixture  to  that  of  a known  gas.  Since, 
for  a stationary,  vertical  column  of  fluid  of  height  h,  the  pressure 
at  the  top  of  the  column  (p j ) is  related  to  the  pressure  at  the 
bottom  (p^ ) by  the  equation 

p2  = Pj  + p gh  , C.  1 

this  comparison  of  hydrostatic  pressures  can  be  made  by  measuring 
the  pressure  difference  at  the  tops  of  two  columns  of  fluid,  one 
being  the  mixture  and  the  other  a known  gas,  and  allowing  the 
pressures  at  the  bottoms  to  be  equal.  Thus,  for  the  mixture 
column, 

P-.  = p1rv,  + P gh 

dm  1 m m 

and  for  the  known  gas  column 

P2g  * »lg  + V • 


Since  ant*  ^2g  are  arran8e<*  to  he  equal. 


The  test  column  can  be  calibrated  in  the  sense  of  determining  gh 
by  comparing  two  known  gases  in  the  above  manner.  Thus, 


where  the  subscript  c is  used  to  denote  calibration  conditions 


Combining  Eqs.  C.  2 and  C.  3 gives 


which,  upon  solving  for  p becomes 


In  the  actual  measurements,  the  gas  that  the  mixture  was  compared 
with  was  air  and  the  calibration  gas  was  helium.  The  pressures 
at  the  bottoms  of  the  columns  were  arranged  to  be  equal  simply 
by  opening  the  bottoms  to  the  atmosphere.  This  was  a stable  con- 
figuration since  the  mixture  was  lighter  than  air.  A column  was 
needed  for  the  air  to  eliminate  room  currents.  In  order  to  insure 


at  least  0.  1%  accuracy  in  the  measure  of  density,  the  temperatures 
of  the  gases  and  the  humidity  of  the  air  had  to  be  monitored  during 
both  the  calibration  and  test  procedures.  Taking  this  into  account, 
Eq.  C.  4 becomes 


where  subscript  T signifies  test  conditions. 

This  procedure  was  suggested  by  Dr.  John  LaRue  of  the 
University  of  California,  San  Diego. 
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APPENDIX  D 

Analysis  for  Predicting  Chemical  Reaction  Rates  from  Scalar 
Probability  Density  Functions 

The  idea  of  using  the  probability  density  functions  of  scalar 

contaminates  in  nonreacting  flows  to  predict  reaction  rates  for 

2 3 

similar  reacting  flows  was  first  proposed  by  Toor  . The  tech- 
nique of  applying  these  ideas  to  the  present  work  are  presented 
below.  Some  of  the  ideas  put  forth  in  Reference  24  are  also 
incorporated. 

For  a nonreacting  shear  layer,  it  is  assumed  that  the 
concentration  histories  of  one  of  the  gases  (C^j)  is  known  as  a 
function  of  both  position  and  time.  Therefore,  the  P.  D.  F.  's  of 
the  concentration  (P  [ C__  3)  can  be  computed.  For  a similar  shear 
layer  in  which  reactant  A is  diluted  in  gas  1 and  reactant  B is 
diluted  in  gas  2 (see  Fig.  1)  and  they  react  irreversibly  according 
to 

NAA  + NBB  = NPP*  D*  1 


This  will  eliminate  the  necessity  of  carrying  the  N.'s  along  in  the 
analysis  but  the  results  will  still  be  valid  for  all  values  of  N.. 


Next,  let 


If  the  diffusivities  of  all  species  are  equal  (D^  = = Dp  = Dx  = D3) 

and  if  the  reaction  rate  is  much  father  than  the  diffusion  rate,  then 


This  is  true  since  CA/C 


A similar  analysis  applies  to 
Combining  D.  2 and  D.  3 with 


the  sum  of  C.  and  C—  will  be  C 


Z2  except  that  Cb/CBq  = 1 
D.  4 and  D.  5 gives 


C„,  then  their  P.  D.  F. 's  will  be  equal 


Now,  since  Y = Z 


With  the  assumed  very  fast  chemistry,  A and  B cannot  exist  at 
the  same  place  at  the  same  time,  or 


Therefore,  if 


then  obviously, 

CA  = X if  X > 0 
= 0 if  X < 0 


Thus, 


and 


cB  = - X if  X < 0 
= 0 if  X > 0 . 

PA[CA]  = PXCCA3  for  CA  > 0 
PB[CB]=  ^ - CB]f°r  CB>0  ' 


From  the  identity 


PX[X]  dX  = PY[Y(X)]dx||  , 


it  can  be  shown  that 


P Jr, 


X + C 


Bo 


YL  C » + C 

■ i.A°t  c-p- 


Ao  Bo 


using  8Y/8X  = 1/(CAq  + CBq),  D.  6 and  D.  9. 
Therefore, 


pa[ca]  = p„LCa  i ^"q1/(C^  + Cnj 


' ' jt. 


.-V  ■ >. 


/<CAo+CBo>‘  D-16 


Since  P„  is  known,  P.  and  P can  be  found.  Thus,  using  the 

il  A D 


identities 


C.  s / c.p.rc.i  dC.  , 

1 J 1 iL  i 1 i 

C. 

1 

c2  ~ I c.sP.rc.i  dc.  , 

i J i lJ  1 


D.  17 


D.  18 


C.'2  = C?  - (&.)  , 

i i r 


D.  19 


CA'  CB* 


and  can  be  calculated. 


Equation  D.  1 5 gives  the  P.  D.  F.  for  for  all  time  that 

C.  >0.  There  is  a delta  function  at  C.  =0  which  accounts  for 
A A 

all  the  time  that  B is  in  excess  making  CA  = 0.  Although  the  area 
of  this  delta  function  can  easily  be  calculated,  it  is  not  necessary 
to  do  so  since  the  P.  D.  F. ’s  are  only  used  as  integrands  when 
multiplied  by  their  arguments  (see  D.  1 7 and  D.  18).  Thus, 

CAP[CA]  = 0 for  C.  = 0.  As  an  example  of  the  use  of  D.  15, 
consider  the  following  given  shape  for  an(*  let  ^Ao  = ^Bo  = 

0.  5.  Then 


C +C 

PAtCAl  = C— nr  - PH  CTT^1  ° < CA  5 CAo 

Ao  Bo  Ao  BoJ 


= PhCCa  + 0.  5]  for  0 < CA  s 0.  5 . 
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Delta  Function 


x Area  Of  Delta 
.Hj  Function*  ^ 


In  order  to  find  C^,  D.  2 and  D4  must  be  used.  Thus 
CA  + CP  . 


From  D.  19, 


c = cA  c„  - CA  . 

p Ao  H A 


= Cp2  - (Cp)2 


Thus,  in  order  to  find  Cp,a,  Cp8  must  be  calculated.  From  D.  20 


Cp3  = CAo2^-  2CAo  Z^CA  + CA3  • 


Since  Zj3  = Cp3  and  CA2  are  known,  ZxCA  must  be  determined. 
From  P.  D.  F.  analysis, 

^"a  = / Zl  (CA>  CA  PA  [CAldCA  • 
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i'i 

In 

It; 


then, 


Z,C 


Z'  'CA>  ° t?  t C°°  for  CA>0- 
Ao  Bo 


CA  + C.  C A + C_ 

C / A BO  w A \n  r A Bo  -i 

J (c;~ ; c~')(cr  ~ cp  )ph[ca  + c " ^dcA» 

„ Ao  Bo  Ao  Bo  Ao  Bo 


D.  24 


or, 


Z’CA  = C-^C 
Ao 


CB»rA 


Bo  CAo  + CBo 


D.  25 


D.  26 


Therefore, 


C°  ■ <=aI  ^ 


D.  27 


C.  CD  CT  2 

Ao  Bo  A /— 

" c + c_  _'cp' 

Ao  Bo 

It  is  interesting  to  note  that  the  integral  definition  of 
unmixedness  takes  on  physical  significance  when  viewed  in  light  of 
this  analysis.  This  integral  definition  is,  restated, 

./(C-  C)dt1+  /(C-  Qdts 
UM  = k 


D.  28 


/<*  - C)dtx  + J C dta 


where  ti  is  the  time  for  C > ^ and  is  the  time  for  C < C. 

From  the  definition  of  probability  densities,  it  can  be  seen  that 

P^j[C]  dC  = dt/T 

where  T is  the  total  time  of  the  sample.  This  simply  states  that 
Pj^tC]  dC  is  the  fraction  of  time  that  C is  within  the  limits  of  dC. 


'•-A. 


\~y 

r 


a 

1 
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Therefore,  D.  28  can  be  rewritten  as 

_/  (C-C)Pjj[C]  dC  + j (C-OPj^C]  dC 


UM  = 


Pj^CldC 


+ c 

cJ  0 


f 


D.  29 


ph[C] 


dC 


These  integrals  can  be  evaluated  by  the  following  analysis.  If  the 
freestream  concentrations  are  arranged  for  each  C such  that 

C, 


Bo 


C . + C _ 

A o Bo 


= C 


D.  30 


then,  whenever  C = C,  no  reactants  will  be  present,  only  product. 

The  reactants  will  have  exactly  "burned"  each  other  up.  If  C > C, 

then  A will  be  in  excess  and  C - C = C . /(C . +C„  ).  If  C < C. 

A ' Ao  Bo'  ' 

then  B will  be  in  excess  and  C - C = CD/(C„  +C„  ).  Using  these 

B Ao  Bo  e 

two  relationships,  along  with  D.  30,  D.  29  becomes 


UM 


0 


r Ao 

Ca  p 

Ca  + Cn 

A Bo 

CAo+CBo  H 

(f 

Ca  +CD 

A o Bo 

dC 


"'Ao  JPH^dC  +CBo{  PHCC]dC 


r Bo 

Cb  r 

"CB+CBo' 

dCB 

CA  +C„  H 

Ao  Bo 

w 

CAo+CBo 

r 1 

& 

CAo  1 prfC]dc  + CBo 

C 

0 

From  Eqs.  D.  15  and  D.  16,  it  can  be  seen  that  the  numerator  is 
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simply  the  sum  of  the  concentrations  of  the  unburned  reactants 


(CA  + C„)  and  the  denominator  is  simply  the  freestream  concentra- 


tions  multiplied  by  the  percentages  of  time  that  each  would  have 
occurred  if  the  gases  were  perfectly  unmixed.  Therefore,  UM 
finally  can  be  interpreted  as 


UM  = 


CA  + CB 


D.  31 


^CA*UM  +(CB)UM 


where  (C.)TTW  is  the  mean  concentration  that  would  have  been 
l UM 


measured  if  the  gases  were  perfectly  unmixed.  Thus,  if  the  gases 
are  perfectly  mixed  (C(t)  = C for  all  time),  then  = Cg  = 0 and 
UM  = 0.  Further,  if  the  gases  are  perfectly  unmixed  (C(t)  = 1 

and  UM  = 1 . 


or  0),  then  C ^ = (CA)UM  and  CB  = (Cn) 


B'UM 


? 


l 


I 


. 

\ 

I 

► 

f#» 
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